US  Army  Corps 
of  Engineers® 

Engineer  Research  and 
Development  Center 


Assessment  of  Effects  of  Maneuver  Training 
Activities  on  Red-cockaded  Woodpecker 
Populations  on  Fort  Stewart,  GA 

Timothy  J.  Hayden,  Robert  H.  Melton,  Beth  Willis,  August  2002 

L.B.  Martin  III,  and  Tim  Beaty 


□  1999 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE  Form  Approved  omb  No. 

0704-0188 

Public  reporting  burder  for  this  collection  of  information  is  estibated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing 
and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burder  to  Department  of  Defense,  Washington 
Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of 
law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (FROM  -  TO) 

01-08-2002  Final  xx-xx-2001  to  xx-xx-2002 

4.  TITLE  AND  SUBTITLE 

Assessment  of  Effects  of  Maneuver  Training  Activities  on  Red-cockaded  Woodpecker 
Populations  on  Fort  Stewart,  GA 

Unclassified 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Hayden,  Timothy  J.  ;  Author 

Melton,  Robert  H.  ;  Author 

Willis,  Beth  ;  Author 

Martin,  L.B.  Ill ;  Author 

Beaty,  Tim  ; 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.S.  Army  Engineer  Research  and  Development  Center  (ERDC) 

PO  Box  9005 

Champaign,  IL6 1826-9005 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

ERDC/CERL  TR-02-17 

9.  SPONSORING/MONITORING  AGENCY  NAME  AND  ADDRESS 

5 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

APUBLIC  RELEASE 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

In  RCWs,  several  factors  have  been  identified  that  affect  reproduction  and  cluster  occupancy  include  both  demographic  and  habitat 
parameters. 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 
Same  as  Report 
(SAR) 

18. 

NUMBER 
OF  PAGES 
74 

19.  NAME  OF  RESPONSIBLE  PERSON 

Hayden,  Timothy 

Timothy.J.Hayden@erdc.usace.army.mil 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

Unclassified  Unclassified  Unclassified 

19b.  TELEPHONE  NUMBER 

International  Area  Code 

Area  Code  Telephone  Number 

217352-6511 

DSN 

Standard  Form  298  (Rev.  8-98) 
Prescribed  by  ANSI  Std  Z39.18 


2 


ERDC/CERL  TR-02-17 


Foreword 


This  study  was  conducted  for  Headquarters,  Department  of  the  Army,  Assistant 
Chief  of  Staff  (Installations  Management)  (ACS(IM))  under  62720A896,  "Envi¬ 
ronmental  Quality  Technology,"  Work  UnitTDO,  'Threshold  I  mpact  of  M aneuver 
Training.”  Thetechnical  monitor  was  Dr.  Victor  Diersing,  DAIM-ED-N. 

The  work  was  performed  by  the  Ecological  Processes  Branch  CN-N  of  the  I  nstal- 
lation  Division  (CN),  Construction  Engineering  Research  Laboratory  (CERL). 
TheCERL  Principal  Investigator  wasTimothyJ  .  Hayden.  Robert  H.  Melton  is  a 
researcher  employed  by  CERL.  Beth  Willis  and  Tim  Beaty  are  employed  by  the 
DPW,  Fish  and  Wildlife  Branch  at  Fort  Stewart,  GA.  L.B.  Martin  III  is  em¬ 
ployed  by  the  Department  of  Ecology  and  Evolutionary  Biology  at  Princeton 
University,  Princeton,  NJ  .  Thetechnical  editor  was  Gloria  J  .  Wienke,  Informa¬ 
tion  Technology  Laboratory.  Steve  Hodapp  is  Chief,  CEERD-CN-N,  and  Dr. 
J  ohn  T.  Bandy  is  Chief,  CEE RD-CN.  The  associated  Technical  Director  was  Dr. 
William  D.  Severinghaus,  CEERD-CV-T.  The  Acting  Director  of  CERL  is  Dr. 
Alan  W.  Moore. 

CERL  is  an  element  of  the  U.S.  Army  Engineer  Research  and  Development  Cen¬ 
ter  (ERDC),  U.S.  Army  Corps  of  Engineers.  The  Commander  and  Executive  Di¬ 
rector  of  ERDC  is  COL  J  ohn  W.  Morris  III,  EN  and  the  Director  of  ERDC  is  Dr. 
J  ames  R.  Houston. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional 
purposes.  Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of 
such  commercial  products.  All  product  names  and  trademarks  cited  are  the  property  of  their  respective 
owners.  The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position 
unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  IT  IS  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 


ERDC/CERL  TR-02-17 


3 


Contents 

Foreword . 2 

List  of  Figures  and  Tables . 5 

1  Introduction . 7 

Background . 7 

Objective . 8 

Approach . 8 

Scope . 10 

Mode  of  Technology  Transfer . 1 0 

2  Methods . 11 

Null  Hypotheses . 11 

Sample  Cluster  Selection . 11 

Training  Data  Collection . 1 2 

RCW  Demographic  Data . 17 

Habitat/Site  Characterization  Data . 17 

Statistical  Analyses . 1 8 

3  Summary  of  Observed  Activity . 1 9 

4  Relationship  of  Range  Scheduling  Data  and  Observations . 25 

Analyses . 25 

Discussion . 28 

5  Data  Reduction . 29 

Demography . 29 

Habitat . 34 

Human  Disturbance . 40 

6  Analyses  of  Maneuver  Training  Effects . 45 

7  Determining  Population  Level  Effects  of  Military  Disturbance . 49 

Materials  and  Methods . 49 

Results . 53 

Discussion . 60 


4 


ERDC/CERL  TR-02-17 


8  Conclusions  and  Recommendations . 62 

Fort  Stewart  Results . 62 

Use  in  Endangered  Species  Risk  Assessments . 64 

References . 67 

Appendix:  PVAvES  Input  Data  File . 69 

CERL  Distribution . 73 

REPORT  DOCUMENTATION  PAGE . 74 


ERDC/CERL  TR-02-17 


5 


List  of  Figures  and  Tables 

Figures 

1  Location  of  sample  clusters . 1 2 

2  Relative  frequency  of  vehicle  type  observed  during  1 0-minute  sample 

observations  during  1997-1999* . 20 

3  Mean  proportion  of  10-minute  sample  observations  in  which  military  activity 

was  observed  during  1997-1999 . 21 

4  Proportion  of  observed  total  human  activity  during  10-minute  sample 

observations  during  1997-1999 . 22 

5  Proportion  of  observed  military  activity  during  10-minute  sample  observations 

during  1997-1999 . 22 

6  Proportion  of  observed  total  human  activity  during  10-minute  sample 

observations  during  1997-99 . 23 

7  Histograms  for  analysis  of  1998  data . 27 

8  Mean  number  of  fledged  by  year,  grouped  by  independent  variables . 32 

9  Ordination  of  RCW  clusters  along  first  three  principal  components  of  measured 

habitat  variables . 37 

1 0  Mean  number  of  young  fledged  by  year,  grouped  by  SITELOAD . 42 

1 1  Mean  number  of  young  fledged  by  year,  grouped  by  ROAD(CLASS) . 43 

12  Classification  of  Red-cockaded  Woodpecker  clusters  used  in  the  population 
viability  analysis,  based  on  2-  or  3-year  mean  levels  of  military  and  civilian 

activity . 54 

13  Estimated  pseudoextinction  probabilities  for  the  Red-cockaded  Woodpecker 

population  on  Fort  Stewart . 58 

1 4  Estimated  probabilities  that  the  Red-cockaded  Woodpecker  population  on  Fort 

Stewart  should  be  classified  as  vulnerable,  endangered,  and  critical . 59 

Tables 

1  Description  and  coding  for  variables  included  for  analyses  or  data  collected  on 

Fort  Stewart,  Georgia  in  1997-1999 . 14 

2  Observed  vehicle  types  during  10-minute  sample  observations,  1997-1999 . 20 

3  Average  percentage  of  10-minute  observation  periods  with  observed  events  per 

monitored  cluster . 20 

4  Events  of  human  activity  observed  in  proximity  to  RCW  clusters . 24 


6 


ERDC/CERL  TR-02-17 


5  Spearman  rank  correlation  matrix  for  military  training  variables  and  cluster 

proximity  to  roads . 26 

6  Non-parametric  means  test  for  NFLDG  grouped  by  demographic  variables . 32 

7  Significance  values  for  Fisher’s  exact  test  of  independence  between  dependent 

variables  (rows)  and  independent  variables  (columns)  1997-1999 . 34 

8  Critical  values  from  principal  component  analysis  of  RCW  cluster  habitat 

variables . 38 

9  Correlation  coefficients  of  PCA  scores  and  measured  RCW  cluster  habitat 

variables  after  varimax  rotation . 38 

1 0  Spearman  rank  correlation  coefficients  between  principal  components  scores 
of  measured  habitat  variables  and  number  of  useable  RCW  cavities  in  1998 

and  number  of  fledglings  produced  by  cluster  for  years  1997,  1998,  and  1999 . 39 

11  Critical  values  (Wald  statistic  (W)  and  P-value  (P))  from  logistic  regression 

analyses  of  principal  component  scores  of  habitat  variables  and  number  of 
useable  RCW  cavities  in  1998  and  binomial  measures  of  RCW  reproductive 
success  in  the  years  1997-1999 . 40 

12  Spearman  rank  correlation  matrix  for  potential  disturbance  variables  and 

number  of  young  fledged,  1997-99 . 41 

13  Non-parametric  means  test  for  NFLDG  grouped  by  SITELOAD  and  ROAD 

(CLASS) . 42 

14  Logistic  regression  statistics . 44 

1 5  Results  of  multiple  linear  regression  for  NFLDG  for  all  independent  variables 

entered  as  model  terms . 46 

1 6  Comparison  of  relative  significance  of  regression  coefficients  of  independent 

variables  in  stepwise  regression  (GLM)  for  number  of  young  fledged . 47 

17  Comparison  between  complete  multivariate  regression  model  (GLM)*  results 

for  NFLDG  with  all  possible  submodels  with  two  independent  variables . 47 

1 8  Nest  productivity  summary  statistics,  over  three  years,  for  Red-cockaded 

Woodpecker  clusters  in  “low”  and  “high”  disturbance  activity  classes . 55 

1 9  Estimated  PVAvES  nest  productivity  parameters,  calculated  from  the  RCW 

clusters  in  “low”  and  “high”  disturbance  activity  classes . 55 

20  Parameterization  of  nine  hypothetical  disturbance  scenarios,  derived  as 

weighted  averages  of  the  parameter  estimates  in  Table  19 . 56 


ERDC/CERL  TR-02-17 


7 


1  Introduction 


Background 

Numerous  U.S.  military  installations  provide  habitat  for  significant  populations 
of  Federally  listed  threatened  and  endangered  species.  I  n  a  1998  survey  of  Army 
installations,  95  installations  reported  on-site  occurrence  of  Federally  listed  spe¬ 
cies  (Schreiber  and  Reed  1999).  The  potential  conflicts  between  compliance  with 
the  Endangered  Species  Act  of  1973,  as  amended  (ESA)  and  the  Army's  mission 
to  train  and  prepare  soldiers  for  military  conflict,  places  a  significant  burden  on 
the  Army  to  evaluate  effects  of  its  mission  activities  on  endangered  species  oc¬ 
curring  on  installation  lands  and  to  support  recovery  of  endangered  populations. 
Efforts  to  resolve  conflicts  between  endangered  species  protection  and  recovery 
and  the  Army’s  training  requirements  have  often  been  hampered  by  a  lack  of 
knowledge  of  the  potential  effects  of  military  training  activities  on  endangered 
species  (Gutzwiller  and  Hayden  1997,  Hayden  1997).  Without  this  information 
the  U.S.  Fish  and  Wildlife  Service  (USFWS),  which  is  the  agency  for  federal 
regulatory  actions  under  the  ESA,  can  elect  to  place  conservative  restrictions  on 
military  training  activities  to  meet  the  mandate  for  protection  and  recovery  of 
endangered  species  populations  and  habitats.  The  negative  impact  of  restric¬ 
tions  on  unit  readiness  due  to  the  presence  of  endangered  species  was  high¬ 
lighted  by  the  Army’s  testimony  before  the  U.S.  Senate's  Environment  and  Pub¬ 
lic  Works  Committee  during  hearings  in  1995  (Hayden  1997). 

As  a  result  of  increasing  conflicts  between  endangered  species  occurrence  and 
training  and  increasing  demand  for  information  on  the  actual  effects  of  training 
on  endangered  species  populations  and  habitats,  the  U.S.  Army  Corps  of  Engi¬ 
neers  (USACE)  initiated  research  in  fiscal  year  1996  under  the  USACE  direct- 
funded  research  program  to  evaluate  effects  of  maneuver  training  on  endangered 
species.  Concurrent  with  the  beginning  of  this  research  initiative,  the  Army  ini¬ 
tiated  a  revision  of  the  Army  management  guidelines  for  the  red-cockaded  wood¬ 
pecker  (Picoides  borealis,  RCW)  on  Army  installations.  This  species  occurs  on 
seven  Army  installations,  including  three  major  Forces  Command  (FORSCOM) 
and  two  major  Training  and  Doctrine  Command  (TRADOC)  installations.  Train¬ 
ing  restrictions  on  these  installations  were  the  impetus  for  revision  of  the  Army 
management  guidelines  to  reduce  restrictions  on  training  while  enhancing  re¬ 
covery  activities  for  the  RCW.  However,  knowledge  of  the  effects  of  maneuver 
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training  on  RCWs  was  lacking,  and  the  USFWS  required  the  Army  to  implement 
a  research  program  to  evaluate  training  effects  on  RCWs.  The  USFWS  docu¬ 
mented  this  requirement  in  its  October  1996  opinion  for  implementing  the  re¬ 
vised  Army  guidelines.  This  requirement  under  the  USFWS  October  biological 
opinion  provided  the  focus  for  the  USACE  direct-funded  research  to  evaluate  ef¬ 
fects  of  maneuver  training  activities  on  the  RCW. 

I  n  early  1997  a  draft  research  plan  was  submitted  to  USFWS  for  peer  review  of  a 
study  to  evaluate  effects  of  military  training  on  RCWs  on  Fort  Stewart,  GA,  that 
would  meet,  in  part,  requirements  for  implementing  the  revised  Army  manage¬ 
ment  guidelines  for  RCWs.  This  plan  (Hayden  1999)  outlined  a  research  ap¬ 
proach  to  evaluate  the  relative  effects  of  military  training,  habitat,  and  demo¬ 
graphic  parameters  on  critical  RCW  population  parameters.  The  report 
documents  the  results  of  these  studies  performed  during  the  years  1997,  1998, 
and  1999. 


Objective 

The  purpose  of  this  research  is  to  evaluate  effects  of  military  maneuver  training 
activities  on  RCW  populations  on  Fort  Stewart,  GA,  and  based  on  these  research 
results  to  provide  recommendations  to  the  Army  for  monitoring  and  assessment 
of  training  effects  on  endangered  species. 


Approach 

Hayden  (1999)  discusses  study  design  considerations  and  limitations  for  this  re¬ 
search  effort.  A  primary  challenge  in  evaluating  potential  impacts  of  maneuver 
training  on  RCW  populations  and  habitats  is  characterizing  and  quantifying 
maneuver  training  in  an  environmental  context.  "Maneuver  training"  is  not  a 
specific  entity  that  can  be  easily  quantified  and  described.  It  is  a  complex  inter¬ 
action  of  events,  participants,  and  equipment  that  is  highly  dynamic  and  vari¬ 
able  both  spatially  and  temporally.  Another  major  challenge  of  any  study  im¬ 
plemented  at  the  landscape  scale  is  quantifying  the  temporal,  spatial,  and 
stochastic  variability  of  the  natural  system  under  consideration. 

The  high  degree  of  dynamic  variability  both  in  the  natural  system  (RCW  popula¬ 
tions  and  habitats)  and  the  factors  (maneuver  training  and  RCW  management) 
potentially  affecting  these  systems  presents  several  difficulties  in  developing  ap¬ 
proaches  to  evaluate  potential  impacts  of  maneuver  training  on  populations  or 
habitats.  First,  the  spatial  and  temporal  variability  of  maneuver  training  as  it  is 


ERDC/CERL  TR-02-17 


9 


conducted  under  actual  conditions  has  not  been  well  documented  and  currently 
cannot  be  easily  predicted  at  the  site-specific  level  —  there  is  no  baseline  data  to 
establish  experimental  levels  for  training  that  reflect  predicted  or  anticipated 
levels  of  training  activity.  Second,  given  the  inherent  complexity  in  the  charac¬ 
teristics  of  maneuver  training,  extreme  care  must  be  given  to  constructing  ap¬ 
propriate  hypotheses  that  can  be  evaluated  given  the  complexity  of  the  system 
and  factors  under  consideration.  Third,  results  of  any  observed  effects  at  the 
site-specific  or  individual  level  must  be  evaluated  in  the  context  of  population 
effects. 

Data  necessary  to  address  these  research  complexities  and  to  test  established 
hypotheses  can  be  derived  from  either  of  two  fundamental  research  approaches: 
experimental  designs  or  observational  studies.  Hayden  (1999)  discusses  several 
factors  that  preclude  controlled  experimental  designs  to  evaluate  hypotheses 
specified  in  this  study.  Limitations  in  implementing  an  experiment  are  due  pri¬ 
marily  to  the  complexity  of  mission  activities  related  to  maneuver  training.  Due 
to  the  factors  discussed  in  Hayden  (1999),  the  current  study  relies  on  observa¬ 
tional  data  to  evaluate  effects  of  training  on  RCW  populations  on  Fort  Stewart. 
Data  collection  and  analyses  emphasize  three  major  areas  of  investigation: 

•  Characterizing  maneuver  training  activities  in  relation  to  RCW  populations 
and  habitats. 

•  Evaluating  effects  of  maneuver  training  activities  on  critical  RCW  population 
parameters  in  relation  toother  habitat  and  demographic  variables  affecting 
RCW  populations. 

•  Evaluating  observed  or  potential  effects  of  training  at  the  population  level. 

This  report  first  documents  research  methods  and  analytical  approaches  in 
Chapter  2.  Chapter  3  provides  data  summaries  for  observations  of  maneuver 
training  activities  in  proximity  to  RCW  clusters  and  Chapter  4  evaluates  mili¬ 
tary  scheduling  data  as  a  predictor  of  site-specific  training  activities.  Chapters  5 
and  6  report  the  results  for  field  studies  conducted  on  Fort  Stewart  during  1997- 
1999  to  assess  effects  of  training  on  RCW  populations.  Chapter  7  demonstrates 
application  of  results  from  field  studies  to  evaluate  extinction  risk  and  recovery 
probability.  Chapter  8  reports  conclusions  from  research  on  Fort  Stewart  and 
discusses  applications  for  generalized  approaches  to  assessment  of  training  ef¬ 
fects  on  endangered  species  and  military  installations. 
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Scope 

This  report  documents  research  results  for  field  studies  conducted  on  Fort  Stew¬ 
art,  GA,  during  1997-1999.  Conclusions  of  this  report  based  on  empirical 
observations  specific  to  data  collected  on  Fort  Stewart  should  not  be  extrapolated 
toother  locations.  General  research  and  analytical  approaches  and  applications 
of  empirical  models  may  be  applied  to  assessment  of  military  training  effects  in 
other  locations  and  contexts. 


Mode  of  Technology  Transfer 

This  report  will  be  made  accessible  through  the  World  Wide  Web  (WWW)  at 
URL: 

http://www.cecer.armv.mil 
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2  Methods 

Based  on  the  research  design  considerations  discussed  in  Chapter  1  and  Hayden 
(1999),  the  approach  for  this  research  relied  primarily  on  observational  data  to 
investigate  the  relationship  of  maneuver  training  activities  and  RCW  demo¬ 
graphic  parameters.  The  null  hypotheses  outlined  here  are  formulated  based  on 
the  availability  of  observational  data  to  test  measures  of  association.  This  ap¬ 
proach  in  general  will  not  establish  cause-and-effect  relationships  between  ma¬ 
neuver  training  activity  and  RCW  demographic  and  habitat  parameters.  How¬ 
ever,  this  approach  is  designed  to  provide  adequate  information  to  make 
informed  evaluations  and  decisions  regarding  the  associations  of  training  activ¬ 
ity  and  RCW  populations  relative  to  effects  attributable  to  habitat  and  demo¬ 
graphic  factors.  This  chapter  revises  proposed  methods  for  this  research  de¬ 
scribed  by  Hayden  (1999). 


Null  Hypotheses 

The  following  generalized  null  hypotheses  (HO)  guided  data  collection,  summary, 
and  statistical  analyses  to  characterize  maneuver  training  activity  in  RCW  clus¬ 
ters  and  evaluate  relationships  between  maneuver  training  activity  and  RCW 
populations  and  habitats. 

•  No  relationship  exists  between  levels  of  training  activity  indicated  from 
range  scheduling  data  and  observations  of  training  activity  at  cluster  sites. 

•  No  relationship  exists  between  RCW  demographic  parameters  and  maneuver 
training  activity  on  Fort  Stewart,  GA,  during  1997-1999. 

These  general  hypotheses  are  further  specified  as  statistical  tests  presented  in 
Chapters  3  through  6. 


Sample  Cluster  Selection 

All  training,  demographic,  and  habitat  data  evaluated  in  this  report  were  col¬ 
lected  at  RCW  clusters  where  demographic  monitoring  of  RCW  populations  was 
conducted  by  Fort  Stewart  biologists  (Figure  1).  Sample  clusters  for  this  report 
are  located  in  each  of  the  six  training  areas  (A,  B,  C,  D,  E,  and  F).  Maneuver 
training  is  conducted  in  all  training  areas  on  Fort  Stewart.  Sample  clusters  in 
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these  areas  represent  an  approximately  25  percent  random  sample  of  all  RCW 
clusters  occurring  in  these  areas  (Fort  Stewart  Endangered  Species  Management 
Plan,  in  draft).  No  sample  clusters  in  this  report  were  located  in  designated  fir¬ 
ing  ranges  or  impact  areas.  Vehicle  maneuver  and  dismounted  personnel  typi¬ 
cally  do  not  train  in  these  areas  except  in  designated  areas  (e.g.,  on  firing  lines). 


□  1999 

Figure  1.  Location  of  sample  clusters. 


Training  Data  Collection 

Characterizing  training  activity  at  cluster  sites  was  done  using  essentially  a 
point  sample  of  training  activity  in  association  with  monitored  RCW  primary 
clusters.  The  protocol  described  here  requires  minimal  interaction  with  training 
units.  The  advantages  of  this  approach  are  that  it  (1)  characterizes  training  ac¬ 
tivity  under  actual  conditions,  (2)  provides  data  that  are  easily  statistically 
summarized  and  analyzed,  and  (3)  requires  minimal  interaction  with  unit  train¬ 
ing  so  that  the  training  "behavior”  is  not  influenced  by  researcher  observation. 

Sample  Unit 

The  unit  of  sampling  at  each  sample  cluster  was  a  10-minute  sample  observa¬ 
tion.  This  time  period  was  selected  as  the  median  of  the  expected  maximum  du¬ 
ration  of  training  activities  (20  minutes)  in  clusters  based  on  information  pro¬ 
vided  in  the  biological  assessment  of  the  1994  Army  guidelines  (Hayden  1997). 


ERDC/CERL  TR-02-17 


13 


Sample  Size 

T raining  data  from  the  10-min  sample  observations  were  analyzed  for  24  sample 
clusters  from  1997.  Data  were  analyzed  for  51  clusters  from  1998,  and  60  clus¬ 
ters  from  1999.  A  minimum  of  twenty  10-minute  observations  for  each  cluster 
were  required  for  a  cluster  to  be  included  in  training  activity  analyses  in  this  re¬ 
port.  These  clusters  were  primary  clusters.  Although  data  were  also  conducted 
at  recruitment  clusters  during  1997-1999,  these  were  not  included  in  compara¬ 
tive  analyses  because  recruitment  clusters  were  at  variable  stages  of  manage¬ 
ment  intervention  (e.g.,  cavity  availability)  and  restricted  status  and  would  not 
be  directly  comparable  to  primary  clusters. 

Sampling  Period 

All  10-minute  sample  observations  were  collected  during  the  period  1  April 
through  30J  une  in  1997-1999.  This  calendar  period  coincides  with  the  peak  sea¬ 
son  for  RCW  reproductive  and  nesting  activities. 

Sampling  Protocol 

The  following  protocol  was  developed  to  balance  sampling  efficiency  with  tempo¬ 
ral  and  spatial  randomization  and  representation  of  sampling  observations. 

•  Sample  clusters  wereallocated  into  sample  groups  to  increase  travel  effi¬ 
ciency  between  cluster  sites. 

•  Observations  in  sample  groups  were  conducted  in  one  of  three  diurnal  peri¬ 
ods:  Morning-  sunrise  to4h  post-sunrise;  Mid-Day-  4h  post-sunrise  to  4  h 
before  sunset;  Evening  -  4  h  before  sunset  to  sunset. 

•  Observations  in  sample  groups  were  rotated  among  diurnal  periods.  For  ex¬ 
ample,  on  Day  1  sample  group  "A"  was  sampled  during  the  "Morning”  period 
and  sample  group  "B"  was  sampled  during  the  "M  id-Day”  period.  On  Day  2, 
sample  group  "A"  would  be  sampled  during  the  "Mid-Day"  period  and  sample 
group  "C"  would  be  sampled  during  the  "Evening"  period,  etc. 

•  Each  sampling  rotation  was  completed  prior  to  the  start  of  a  subsequent 
sampling  rotation  to  ensure  that  all  clusters  were  sampled  an  equal  number 
of  ti  mes. 

Training  Data  Recorded  at  Sample  Sites 

Table  1  provides  detailed  descriptions  and  coding  for  training  variables  recorded 
at  sample  clusters  and  used  for  analyses  in  this  report.  During  each  10-minute 
observation  period,  an  observer  parked  or  walked  to  a  cluster  location,  but  out¬ 
side  the  200-foot  (61-m)  radius  restricted  area  for  military  training  (Hayden 
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1997)  and  recorded  data  on  all  human  activity  in  proximity  to  the  cluster.  Hu¬ 
man  "events"  were  recorded  if  an  activity  was  in  the  visual  perceptual  range  of 
an  RCW  at  the  cluster  site,  which  would  be  inclusive  of  all  areas  within  the  200- 
foot  (61-m)  buffer  and  would  potentially  include  observed  activities  beyond  the 
200-foot  (61-m)  buffer.  Although  this  required  a  somewhat  subjective  judgment 
on  the  part  of  the  observer,  this  determination  was  usually  unambiguous  due  to 
the  nature  of  the  activities  and  the  proximity  of  these  activities  to  clusters.  Data 
recorded  for  events  observed  during  a  10-minute  observation  period  included 
type  of  personnel  (civilian,  civilian  staff,  or  military),  number  and  type  of  vehi¬ 
cles  and  personnel,  duration  of  event,  nature  of  event  (off-road  or  road  transit), 
and  notes  on  any  other  characteristics  of  the  observed  event.  More  than  one 
event  could  occur  during  a  10-minute  observation  period.  An  event  was  defined 
as  any  related  set  of  activities  by  a  group  of  personnel  or  vehicles.  For  example, 
a  convoy  of  six  vehicles  in  road  transit  would  be  recorded  as  one  event.  A  visit  by 
Fort  Stewart  Fish  and  Wildlife  Branch  personnel  observed  during  the  same  10- 
minute  period  would  be  recorded  as  a  separate  event.  Data  recorded  included 
type  of  activity  observed,  number  of  troops  and/or  vehicles  involved,  duration  of 
the  activity,  and  personnel  type  (military,  civilian,  or  civilian  staff  of  Fort  Stew¬ 
art). 


Table  1.  Description  and  coding  for  variables  included  for  analyses  or  data  collected  on  Fort 
Stewart,  Georgia  in  1997-1999. 


Variable  name 

Description 

Coding 

Dependent 

variables 

NFLGD 

Number  of  young  fledged  /  breeding  pair 

SUCCESS 

Were  young  fledged  at  the  cluster 

0  =  No 

1  =  Yes 

BREEDING 

Was  a  breeding  pair  present 

0  =  No 

1  =  Yes 

MRET 

Breeding  male  retention. 

Y  =  breeding  male  same  as  previous  year 

N  =  breeding  male  from  previous  year  not 
present 

FRET 

Breeding  female  retention 

Y  =  breeding  female  same  as  previous  year 

N  =  breeding  female  from  previous  year  not 
present 

MREC 

Breeding  male  recruitment 

Y  =  Recruitment  of  new  breeding  male 

N  =  No  recruitment  of  breeding  male 

FREC 

Breeding  female  recruitment 

Y  =  Recruitment  of  new  breeding  female 

N  =  No  recruitment  of  breeding  female 
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Variable  name 

Description 

Coding 

Demographic 

variables 

GRPN 

Presence  of  auxiliaries  with  breeding  pairs 

Y  =  auxiliaries  are  present  in  group  with 
breeding  pairs 

N  =  auxiliaries  not  present 

PRSUC 

Whether  young  were  fledged  at  the  cluster 
in  the  previous  year 

Y  =  >  0  young  were  fledged 

N  =  0  young  were  fledged 

ISOL 

Two-level  k-means  classification  of  number 

NO  =  lower  level  of  isolation  from  neighbor- 

of  active  clusters  within  two  kilometers  of 
sample  cluster 

ing  clusters 

YES  =  higher  level  of  isolation  from 
neighboring  clusters 

MBRYR 

Year  as  breeding  male  in  a  sampled  cluster 

1  =  first  year  as  breeding  male  in  cluster 
(may  not  be  first  year  as  breeder) 

2  =  second  year  or  more  as  breeding  male 
incluster 

FBRYR 

Year  as  breeding  female  in  a  sampled  clus¬ 

1  =  first  year  as  breeding  female  in  cluster 

ter 

(may  not  be  first  year  as  breeder) 

2  =  second  year  or  more  as  breeding  female 
incluster 

Habitat 

variables 

PINEBA 

Total  pine  basal  area  /  ha  for  cluster 

PINEL10 

Pine  BA/ha  for  stems  <  10"  DBH 

HDBA 

Hardwood  BA/ha  for  cluster 

PINESTM 

pine  stems  /  ha  by  cluster 

HDSTM 

Number  of  hardwood  stems  /ha 

5BURN 

Total  burns  in  5  years  prior  to  current  breed¬ 
ing  season 

IBURN 

Interval  of  most  recent  burn  to  current 
breeding  season 

1=  within  one  year 

2=  1  -2  years 

3  =  2-3  years 

4  =  3-4  years 

5  =  4-5  years 

6  =  >  5  years 
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Variable  name 

Description  Coding 

Training 

variables 

NTRP 

Average  RFMSS  troops/day  scheduled  for 

TA,  Apr-June  1998. 

NDAY 

total  RFMSS  days  scheduled  for  TA,  Apr- 
June  1998. 

ROAD 

Distance  in  meters  from  nest  tree  or  ran¬ 
domly  selected  cavity  tree  to  nearest  main¬ 
tained  road. 

LIOROAD 

Log  transformation  of  ROAD  variable 

ROADCL 

k-means  classification  of  LIOROAD  variable 

ACT 

Proportion  of  10-min  sample  observations  in 
which  human  activity  (military  or  civilian) 
were  observed 

MILACT 

Proportion  of  10-min  sample  observations  in 
which  military  activity  were  observed. 

CIVACT 

Proportion  of  10-min  sample  observations  in 
which  civilian  activity  were  observed. 

SITELOAD 

k-means  classification  of  arcsine  trans¬ 
formed  MILACT  variable 

Descriptor 

variables 

CLUSTER 

Ft.  Stewart  cluster  identification  number 

TA 

Sub  -  training  area  designation 

AREA 

Major  training  area  designation 

YRACT 

Year  the  cluster  was  activated  for  occupa¬ 
tion 

EVTPER 

total  number  of  events  with  dismounted 

personnel 

NUNIT8 

total  RFMSS  units  scheduled  for  TA,  Apr- 
June  1998.  ("UNITNAME"  by  "TRA  AREA 

Training  Data  From  G3  Range  Scheduling  Records 

Scheduling  data  from  the  Army's  Range  and  Facilities  Management  and 
Scheduling  System  (RFMSS)  for  training  ranges  in  which  sample  clusters  were 
located  on  Fort  Stewart  was  used  to  correlate  scheduled  training  activity  with 
field  observations  of  training  activity  at  sample  clusters.  Based  on 
recommendations  of  FORSCOM  and  installation  trainers  (Army  and  USFWS 
meeting,  Atlanta,  GA,  April  2000),  variables  considered  in  this  report  were: 
number  of  scheduled  training  days  by  training  area,  average  number  of  troops 
per  day  scheduled  by  training  area,  and  number  of  units  scheduled  per  training 
area.  RFMSS  data  were  analyzed  for  the  period  1  April  through  30  J  une  1998. 
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Table  1  provides  description  and  coding  for  training  data  variables  derived  from 
RFMSS  for  analysis  in  this  report.  RFMSS  data  for  1997  were  not  directly 
comparable  to  1998  data.  RFMSS  data  for  1999  were  not  available  at  the  time 
this  report  was  written. 


RCW  Demographic  Data 

RCW  demographic  data  for  monitored  sample  clusters  were  collected  by  person¬ 
nel  of  the  Fort  Stewart  Fish  and  Wildlife  Branch,  DPW,  in  accordance  with 
methods  described  in  the  Fort  Stewart  ESMP  approved  by  the  installation  Com¬ 
mander  in  2001.  These  data  have  been  collected  in  monitored  clusters  since 
1994.  Demographic  variables  obtained  from  the  Fort  Stewart  monitoring  pro¬ 
gram  are  described  and  coded  in  Table  1. 


Habitat/Site  Characterization  Data 

Vegetation  and  site  data  were  collected  at  sample  plots  during  J  uly  in  1998  and 
during  J  uly,  August,  and  September  in  2000.  Variable  description  and  coding 
for  habitat  data  are  shown  in  Table  1.  Vegetation  data  collected  at  each  sample 
plot  included: 

•  Stand  data:  Data  were  collected  using  the  point-centered  quarter  method 
(Cottam  and  Curtis,  1956),  at  five  sample  points  within  each  cluster.  Sample 
points  were  centered  at  cavity  or  start  trees  in  most  cases.  I  n  those  cases 
where  non-cavity  or  non-start  trees  were  used,  only  pines  that  fit  the  pa¬ 
rameters  for  i  nstal  lation  of  i  nserts  or  Copeyon  starts  were  chosen.  Basal 
area  of  pines  and  hardwoods  was  recorded  using  a  10X  wedge  prism.  I  n  each 
quadrant  around  the  tree,  the  species,  distance  (to  nearest  0.1  m),  and  dbh 
(to  nearest  0.1  inch)  of  the  nearest  canopy  pine  (over  15  m  in  height),  the 
nearest  midstory  pine  (2  m  to  15  m  in  height),  and  the  nearest  hardwood  (at 
least  2  m  in  height)  all  within  at  least  50  meters  of  the  central  tree  were  re¬ 
corded. 

•  Ground  and  canopy  cover:  Data  were  collected  along  100-m  transects  posi¬ 
tioned  such  that  the  50-m  point  intersected  an  imaginary  perpendicular  line 
5-m  from  each  tree  sampled  above.  Canopy  coverage  and  ground  cover  were 
recorded  every  5  m  along  the  transect  beginning  at  the  5-m  mark.  Canopy 
cover  included  longleaf  pine,  loblolly  pine,  slash  pine,  pond  pine,  hardwood, 
and  no  canopy  detected.  Ground  cover  was  recorded  as  bare  ground,  litter, 
graminoid  (except  wi regrass),  wi regrass,  forb,  duff,  dead  wood,  woody  debris, 
or  road. 
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•  Burn  history:  Burn  history  was  obtained  for  all  sample  clusters  based  on 
Fort  Stewart  records  and  maps. 

•  Cavity  number:  Number  of  natural  and  artificial  cavities  for  all  sampledus- 
ters  was  obtained  from  Fort  Stewart  records. 


Statistical  Analyses 

A  variety  of  software  packages  were  used  to  perform  statistical  analyses  and 
tests  reported  in  this  document  including  SYSTAT,  SPSS,  and  SAS.  Details  of 
the  analyses  are  provided  in  the  chapters  reporting  results  (Chapters  4  through 
7).  A  p-value  less  than  or  equal  to  0.05  is  considered  significant  for  all  test  sta¬ 
tistics. 


ERDC/CERL  TR-02-17 


19 


3  Summary  of  Observed  Activity 

As  described  in  Chapter  2,  10-minute  sample  observations  were  conducted  at 
RCW  sample  clusters  during  1  April  -  30  J  une,  1997-1999  to  evaluate  human 
activity  in  proximity  to  RCW  clusters.  The  number  of  monitored  clusters  with 
more  than  twenty  10-minute  sample  observations  was  24,  51,  and  60  in  1997, 
1998  and  1999,  respectively.  The  increase  to  51  in  1998  was  due  to  including 
monitored  clusters  on  the  east  side  of  Fort  Stewart,  for  which  no  10-minute 
observation  samples  were  obtained  in  1997.  The  increase  to  60  clusters  in  1999 
was  through  the  addition  of  clusters  that  were  first  activated  by  artificial  nest 
cavities  in  1998  and  1999.  The  average  number  of  10-minute  observations  per 
cluster  was  26.1,  37.8,  and  32.6  in  1997, 1998  and  1999,  respectively 

The  relative  frequency  of  vehicle  types  observed  in  1997-1999  are  show  in  Figure 
2.  Total  numbers  and  breakdown  by  vehicle  type  are  shown  in  Table  2.  Note 
that  civilian  and  civilian  staff  vehicles  were  typically  4X4  trucks,  with  the  excep¬ 
tion  of  a  limited  number  of  earth-movers  performing  road  maintenance.  "Light 
wheeled  military"  vehicles  include  Hu mvees  and  4X4  light  utility  trucks.  "Heavy 
wheeled  military"  includes  all  wheeled  vehicles  with  three  or  more  axles. 
'Tracked  military"  includes  all  tracked  vehicle  types  observed  including  Ml 
tanks,  Bradley  AFVs,  M113  APCs,  and  MLRS.  Tracked  military  vehicles  com¬ 
prised  no  more  that  10.7  percent  (1997)  of  the  total  number  of  observed  vehicles 
in  any  year.  The  percentage  of  civilian  and  civilian  staff  vehicles  combined 
ranged  from  22.7  percent  in  1997  to  30.0  percent  in  1998. 

The  proportion  of  human  activity  observed  in  proximity  to  clusters  during  the 
10-minute  observation  periods  was  relatively  low  during  the  years  1997-1999. 
The  mean  percentage  of  observations  in  which  human  activity  was  observed  in 
proximity  to  clusters  never  exceeded  8.5  percent  of  the  total  observation  periods 
(Table  3).  Observations  of  military  activity  when  averaged  across  clusters  by 
year  ranged  from  3.2  percent  in  1998  to  6.0  percent  in  1997  (Table  3,  Figure  3). 
The  proportion  of  observations  of  civilian  activity  in  proximity  to  clusters  was 
similar  to  observations  of  military  activity  (Table  3).  Observations  of  civilian  ac¬ 
tivity  ranged  from  an  average  of  2.8  percent  in  1997  to  5.5  percent  in  1998. 
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Frequency  of  Observed  Vehicle  Types 


Civilian  Civilian  staff  Light  wheeled  Heavy  wheeled  Tracked  military 

military  military 

Vehicle  Types 


□  F97  ■  F98  O  F99 


Figure  2.  Relative  frequency  of  vehicle  type  observed  during  10-minute  sample  observations 
during  1997-1999*. 


Table  2.  Observed  vehicle  types  during  10-minute  sample  observations,  1997-1999. 


Vehicle  Type* 

Number 

1997 

% 

Number 

1998 

% 

Number 

1999 

% 

Civilian 

4 

4.8 

18 

10.2 

10 

3.2 

Civilian  staff 

15 

17.9 

53 

29.9 

62 

20.1 

Light  wheeled  military 

36 

42.9 

62 

35.0 

176 

57.1 

Heavy  wheeled  military 

20 

23.8 

26 

14.7 

40 

13.0 

Tracked  military 

9 

10.7 

18 

10.2 

20 

6.5 

Total 

84 

177 

308 

*  “Civilian”  vehicles  are  nonmilitary  and  nonstaff  vehicles.  “Civilian  staff”  are  vehicles  operated  by  Fort 
Stewart  civilian  personnel.  “Light  wheeled  military”  are  all  two-axled  wheeled  including  Humvees  and 
light-duty  trucks.  “Heavy  wheeled  military”  are  all  military  wheeled  vehicles  with  three  or  more  axles. 
“Tracked  military”  are  all  military  tracked  vehicles  including  Ml  series  tanks,  Bradley  AFVs,  Ml  13  APCs, 
and  MRLS.  Percent  column  is  percent  of  total  number  of  vehicles  for  each  vehicle  type. 


Table  3.  Average  percentage  of  10-minute  observation  periods  with  observed  events  per 
monitored  cluster. 


Mean  percentage  of  observations  with  observed  activity 

Activity  Type 

%  1997 

N* 

SD 

%  1 998 

N 

SD 

%  1 999 

N 

SD 

All  activity 

8.5 

24 

13.0 

6.1 

51 

7.9 

6.2 

60 

10.8 

Military  activity 

6.0 

24 

10.4 

3.2 

51 

5.5 

3.7 

60 

8.7 

Civilian  Activity 

2.8 

24 

4.8 

5.5 

51 

4.2 

3.0 

60 

4.1 

*  N  =  number  of  clusters  monitored.  SD  =  standard  deviation  of  the  mean. 
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Mean  frequency  of  military  observations  by  year 


Year 

Figure  3.  Mean  proportion  of  10-minute  sample  observations 
in  which  military  activity  was  observed  during  1997-1999. 

The  bar  charts  in  Figures  4  through  6  show  the  proportion  of  10-minute  observa¬ 
tions  with  observed  activity  by  cluster  for  total  activity,  military  activity,  and 
civilian  activity,  respectively.  The  clusters  in  these  figures  are  rank-ordered  by 
their  3-year  average  of  the  proportion  of  observations  with  activity.  Note  that 
the  level  of  activity  in  relation  to  individual  clusters  is  relatively  consistent 
across  years;  i.e.,  clusters  with  low  frequencies  of  observed  activity  tend  to  be  low 
for  all  years  and  clusters  with  high  frequencies  tend  to  have  high  levels  of  activ¬ 
ity  in  all  3  years.  The  range  of  frequencies  of  10-minute  observation  periods  with 
observed  total  activity  in  proximity  to  clusters  was  0.00  -  0.53.  The  range  of  ob¬ 
served  frequencies  for  military  and  civilian  activity  was  0.00  -  0.46  and  0.00  - 
0.19,  respectively. 
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Frequency  of  observations  with  civilian  and/or  military 

activity 


□  TOTACT7  ■  TOTACT8  □  TOTACT9 


Clusters  are  ordered  on  the  x-axis  by  average  proportion  of  observations  of  activity  over  the  3-year  pe¬ 
riod. 

Figure  4.  Proportion  of  observed  total  human  activity  during  10-minute  sample  observations 
during  1997-1999. 


Frequency  of  observations  with  military  activity 


Cluster  identification  number 


□  MILACT7  ■  MILACT8  □  MILACT9 


Clusters  are  ordered  on  the  x-axis  by  average  proportion  of  military  observations  of  activity  over  the  3- 
year  period. 

Figure  5.  Proportion  of  observed  military  activity  during  10-minute  sample  observations  during 
1997-1999. 
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Frequency  of  observations  with  civilian  activity 
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Cluster  identification  number 


□  CIVACT7  ■CIVACT8  nCIVACT9 

Clusters  are  ordered  on  the  x-axis  by  average  proportion  of  civilian  observations  of  activity  over  the  3- 
year  period. 

Figure  6.  Proportion  of  observed  total  human  activity  during  10-minute  sample  observations 
during  1997-99. 


The  type  of  activity  was  characterized  as  "off-road"  and  "road  transit"  for  both 
civilian  and  military  events  (Table  4).  Virtually  all  military  events  in  proximity 
to  RCW  clusters  were  road  transit  by  vehicles.  No  events  with  off-road  military 
activity  were  observed  in  proximity  to  clusters  in  1997  and  1999.  Only  3  of  146 
observed  events  in  1998  were  military  off-road  activity.  These  events  were  vehi¬ 
cles  and/or  personnel  parked  off-road  outside  the  200  foot  RCW  cluster  buffer 
zone.  T able  4  also  shows  the  percentage  of  events  that  involved  dismounted  per¬ 
sonnel  (individuals  on  foot)  versus  strictly  vehicular  traffic.  Overall,  dismounted 
activity  was  relatively  low  for  all  years  ranging  from  13.5  percent  in  1997  to  14.4 
percent  in  1998. 

RCW  management  and  monitoring  checks  by  Fort  Stewart  DPW  and  Fish  and 
Wildlife  Branch  personnel  comprised  virtually  all  of  the  off-road  events  observed 
during  the  3  years  of  this  study.  In  1998,  23  of  146  observed  events  (15.8  per¬ 
cent)  were  civilian  off-road  activities  within,  or  in  proximity  to,  RCW  clusters.  In 
1997  and  1999  civilian  off-road  events  comprised  13.5  percent  and  15.5  percent, 
respectively,  of  the  total  observed  events.  These  personnel  were  primarily  per¬ 
forming  duties  associated  with  RCW  monitoring,  management,  and  research  ac¬ 
tivities,  which  was  consistent  with  training  regulations  in  effect  during  the  pe¬ 
riod  of  this  study. 
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Table  4.  Events  of  human  activity  observed  in  proximity  to  RCW  clusters. 


Year 

Total 

Events 

Observed* 

Off-road 

Road-transit 

Dismounted 

Personnel 

Vehicle  Traffic 

1997 

74 

10  (13.5) 

64  (86.5) 

10  (13.5) 

64  (86.5) 

1998 

146 

26  (17.8) 

120  (82.2) 

21  (14.4) 

125  (85.6) 

1999 

168 

26  (15.5) 

142  (84.5) 

24  (14.3) 

144  (85.7) 

*  Total  events  observed”  are  all  activity  events  observed  by  year.  More  than  one  event  may  be  re¬ 
corded  in  each  10-minute  observation  period.  Events  are  classified  as  “off-road”  or  “road-transit”  and 
“dismounted  personnel”  or  “vehicle  traffic.”  “Off-road”  events  are  observations  of  vehicles  or  personnel 
off  of  maintained  roads  or  trails  in  proximity  to  RCW  clusters.  Road  transit  events  are  observations  of 
vehicles  or  personnel  traveling  on  maintained  roads  or  trails.  Dismounted  personnel  are  events  in 
which  dismounted  personnel  were  observed  in  proximity  to  RCW  clusters.  Percentage  of  total  events  is 
in  parentheses.  Vehicle  traffic  is  any  vehicle  activity  without  associated  dismounted  activity. 
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4  Relationship  of  Range  Scheduling  Data 
and  Observations 


Analyses 

RFMSS  data  were  obtained  for  Fort  Stewart  for  1998  to  evaluate  the  relation¬ 
ship  between  scheduled  training  events  and  observations  of  site-specific  training 
activity  during  10-minute  sample  observations  at  sample  cluster  sites.  RFMSS 
data  from  1  April  -  30  J  une  were  selected  coincident  with  the  period  of  collection 
of  10-minute  sample  observations. 

Two  variables  were  derived  from  RFMSS  database  for  this  analysis.  The  first 
was  the  total  number  of  days  scheduled  by  training  area  (NDAY)  during  April- 
J  une.  The  second  was  the  average  number  of  troops  scheduled  per  day  by  train¬ 
ing  area  (NTRP)  during  April-J  une.  These  variables  were  identified  by  Army 
trainers  as  the  most  consistently  available  variables  across  RFMSS  versions,  the 
most  likely  to  be  consistently  entered  in  the  database,  and  the  easiest  to  inter¬ 
pret  relative  to  scheduled  training  loads. 

In  addition,  proximity  to  roads  was  considered  a  potential  factor  in  the  I  i  kel  i  hood 
of  observing  training  activities  in  proximity  to  RCW  clusters.  Data  were  col¬ 
lected  for  each  sample  cluster  for  straight-line  distance  to  the  nearest  main¬ 
tained  road  (ROAD,  Table  1). 

The  dependent  variable  used  in  these  analyses  is  the  proportion  of  10-minute 
sample  observations  with  observed  military  activity  during  April-J  une  1998. 

The  Spearman  rank  correlation  matrix  for  these  variables  is  shown  in  Table  5. 
Strong  correlations  were  identified  between  average  number  of  troops  scheduled 
per  day  (NTRP)  and  total  training  days  (NDAY)  scheduled  (Spearman  rank  coef¬ 
ficient  =  0.983).  No  other  strong  correlations  among  pairs  of  variables  were  iden¬ 
tified. 

The  variables  NDAY  and  ROAD  were  retained  for  multivariate  analyses  of  pre¬ 
dictive  models  for  site-specific  training  levels.  Of  the  correlated  variables 
NDAY8  (NDAY  in  1998)  and  NTRP8,  the  latter  was  excluded  because  the  values 


26 


ERDC/CERL  TR-02-17 


entered  into  RFMSS  represent  the  troop  numbers  of  the  unit  scheduled  for  train¬ 
ing  but  not  necessarily  the  number  of  troops  actually  training,  while  N  DAY  more 
closely  represents  the  number  of  days  troops  were  training  in  the  scheduled 
training  areas. 


Table  5.  Spearman  rank  correlation  matrix  for  military  training  variables  and  cluster  proximity  to 
roads. 


MILACT* 

ROAD 

NDAY 

NTRP 

MILACT 

1.000 

ROAD 

0.135 

1.000 

NDAY 

0.257 

0.148 

1.000 

NTRP 

0.260 

0.151 

0.983 

1.000 

*  All  scheduling  data  and  cluster  site  data  are  for  year  1998.  Variables  in  this  matrix  are  not  trans¬ 
formed.  MILACT  =  proportion  of  10-minute  observations  in  which  training  activities  were  observed. 
ROAD  =  proximity  of  cluster  (nest  cavity  tree  or  nearest  active  cavity  tree)  to  nearest  maintained  road. 
NDAY  =  number  days  with  scheduled  training  events  during  April  -  June  1 998.  NTRP  =  average  num¬ 
ber  of  troops  scheduled  per  day. 

Multivariate  regressions  were  performed  to  evaluate  the  prediction  of  site- 
specific  levels  of  training  from  data  on  range  scheduling  and  cluster  proximity  to 
road.  The  dependent  variable  in  these  regressions  was  the  proportion  of  10- 
minute  sample  observations  in  which  military  activity  was  observed  (MILACT). 
The  independent  (predictor)  variables  in  these  regressions  were  NDAY  and 
ROAD.  Regressions  were  performed  using  the  general  linear  model  (GLM)  pro¬ 
cedure  (Systat,  v9.01).  The  values  for  MILACT  were  arcsine  transformed  and 
ROAD  was  log  transformed  to  improve  normality  of  the  variable  distributions. 
The  following  regression  model  was  evaluated: 

MILACT  =  constant  +  NDAY  (arcsi ne transformed)  +  ROAD  (log transformed) 

The  F-ratio  (3.437)  for  this  regression  was  significant  (p  =0.041)  as  was  the  re¬ 
gression  coefficient  for  NDAY  (t  =2.374,  p  =0.022).  However,  the  multiple  R2 
(0.125)  suggests  this  model  is  a  poor  predictor  for  the  dependent  variable 
MILACT8.  Inclusion  of  the  interaction  term  NDAY*ROAD  only  slightly 
strengthens  this  relationship  (multiple  R2  =0.197).  The  weak  correlation  of  the 
model  with  the  dependent  variable,  examination  of  the  residuals  against  pre¬ 
dicted  values,  and  the  non-normal  distribution  of  the  variables,  particularly 
MILACT  (Figure  7),  suggest  that  linear  regression  of  these  predictor  variables 
may  not  be  appropriate. 
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MILACT8 

(A)  proportion  of  military  activity  observed  during  10-minute  sample  periods  (MILACT,  arcsine  trans¬ 
formed) 


(B)  Number  of  scheduled  days  of  training  by  training  area  (NDAY) 


ROAD98 


(C)  Proximity  of  cluster  to  nearest  maintained  road  (ROAD,  log  transformed). 


Figure  7.  Histograms  for  analysis  of  1998  data. 
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Because  MILACT  is  highly  skewed  and  weighted  toward  zero  values,  MILACT 
was  classified  as  a  categorical  variable  (SITE LOAD)  for  values  of  "High"  or 
"Low"  site-specific  training  loads  using  the  k-means  classification  procedure  to 
split  the  cases  into  two  groups.  This  grouping  procedure  classified  40  cases  as 
"Low"  site  training  loads  (mean  value  of  arcsine  transformed  variable  =  0.10) 
and  11  cases  as  "High"  site-specific  training  load  (mean  value  of  arcsine  trans¬ 
formed  variable  =0.69). 

Logistic  regression  was  performed  with  MILACT  reclassified  as  a  categorical  de¬ 
pendent  variable  with  two  levels  and  the  independent  variables  ROAD8  and 
NDAY8.  Neither  of  the  independent  variables  nor  combination  of  variables  are 
significant  predictors  (p  >0.05  for  all  variables)  for  the  dependent  variable  clas¬ 
sified  as  "High"  or  "Low"  site-specific  training  loads.  Log-likelihood  ratios  (indi¬ 
cator  of  model  'fit")  for  all  main  effects  models  and  models  with  interactions 
were  not  significant  (p  >0.05). 


Discussion 

No  strong  correlations  or  predictor  variables  were  found  between  range  schedul¬ 
ing  data  from  RFMSS  and  site-specific  observations  of  training  activity  in  prox¬ 
imity  to  cluster  sites.  This  lack  of  association  was  not  improved  by  including 
proximity  to  maintained  roads  in  the  regression  models.  This  lack  of  association 
is  likely  due  to  the  "resolution"  of  both  the  scheduling  data  used  and  the  low  pro¬ 
portion  of  observed  military  activity  and  the  heavy  bias  toward  cases  with  low  or 
no  observations  of  activity  during  the  10-minute  observation  periods.  The  vari¬ 
ables  analyzed  from  the  range  scheduling  data  provide  only  a  coarse  index  of 
troop  activity  in  training  areas.  The  spatial  resolution  of  RFMSS  is  at  least  an 
order  of  magnitude  greater  than  the  spatial  resolution  of  cluster  specific  observa¬ 
tions  of  training  activity  (100's  of  hectares  per  training  area  versus  10's  of  hec¬ 
tares  or  less  for  cluster  sites).  The  low  proportion  of  observations  of  activity  from 
the  10-minute  sample  observations  is  problematic  for  using  range  scheduling 
data  or  other  factors  such  as  proximity  to  roads  as  predictors  of  training  activity 
at  a  site.  The  data  suggest  that  military  activity  in  relation  to  a  specific  cluster 
site  is  a  relatively  rare  event  with  inherent  statistical  issues  of  sampling  and 
predicting  rare  events.  Further  analysis  of  data  available  from  RFMSS  is  re¬ 
quired  to  evaluate  its  utility  for  classifying  disturbance  at  the  site  level. 
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5  Data  Reduction 

This  chapter  presents  exploratory  analyses  and  statistical  models  designed  to 
reduce  the  "dimensions”  of  the  extensive  data  set  collected  to  evaluate  effects  of 
maneuver  training  on  RCW  populations.  Over  60  measures  of  habitat,  demo¬ 
graphic,  and  training  parameters  initially  were  considered  in  this  study.  The 
purpose  of  the  analyses  in  this  chapter  is  to  reduce  the  number  of  variables  con¬ 
sidered  in  final  multivariate  regression  models  for  assessment  of  maneuver 
training  effects  on  RCW  populations.  These  analyses  seek  to  identify  redundant 
variables,  correlations  among  related  variables,  significant  associations  between 
dependent  and  independent  variables,  adequacy  of  sample  size,  and  appropriate 
variable  transformations.  Retention  of  variables  for  final  analyses  of  training  ef¬ 
fects  is  based  on  significance  of  association  between  dependent  and  independent 
variables,  adequacy  of  sample  size,  and  how  well  they  represent  critical  habitat, 
demographic,  and  training  parameters.  Related  groups  of  independent  variables 
(demographic,  habitat,  and  training  variables)  are  evaluated  independently  to 
aid  interpretation  of  relationships  among  variables  and  measures  of  association 
with  dependent  variables. 

Table  1  (page  14)  provides  descriptions  and  coding  for  all  dependent  variables 
considered  in  these  analyses.  These  variables  reflect  critical  aspects  of  RCW 
demography  and  provide  measures  of  reproductive  success  (number  of  young 
fledged  [NFLDG]  and  nest  success  [SUCCESS])  and  turnover  of  breeding  adults 
in  cluster  sites  (presence  of  breeding  pair  [BREEDING],  retention  [MALERET, 
FEMRET]  and  recruitment  [MALEREC,  FEMRET]  of  breeding  males  and  fe¬ 
males).  I  ndependent  (predictor)  variables  are  those  parameters  known  to  be  or 
potentially  related  to  RCW  reproductive  success  and/or  adult  retention. 


Demography 

This  section  determines  the  relationship  of  demographic  parameters  to  number 
of  young  fledged,  presence  of  breeding  pair,  nest  success,  breeding  male  and 
female  retention,  and  breeding  male  and  female  recruitment.  Independent 
variables  in  these  analyses  areGRPN,  PRSUC,  ISOL,  MBRYR,  and  FBRYR  (see 
Table  1  for  definitions).  MBRYR  and  FBRYR  are  class  variables  with  two  levels 
(1  =  male/female  first  year  of  breeding  in  cluster,  2  =  male/female  second  year  of 
breeding  in  cluster).  The  range  of  original  values  for  GRPN  (number  of 
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individuals  in  the  group)  and  PRSUC  (how  many  young  were  fledged  in  the 
preceding  year)  are  small  (range  0  -  5  for  GRPN,  0  -  3  for  PRSUC),  so  these 
were  transformed  to  class  variables  with  two  levels  for  these  analyses.  I  SOL 
was  transformed  to  a  class  variable  with  two  levels  using  k-means  clustering  to 
simplify  analysis  with  all  independent  variables  as  class  variables  with  two 
levels.  Treatment  of  all  independent  variables  as  class  variables  simplifies 
statistical  treatment  of  the  variables,  simplifies  comparison  of  results  among 
variables,  and  avoids  problems  associated  with  non-normal  distributions  of  the 
original  variables. 

Figure  8(A)  shows  that  the  mean  number  of  young  fledged  by  breeding  pairs 
without  auxiliaries  was  lower  in  all  years  than  breeding  pairs  with  auxiliaries. 
Mean  number  of  young  fledged  was  also  consistently  lower  in  clusters  where 
young  were  not  fledged  in  the  previous  year  versus  clusters  where  young  were 
fledged  in  the  previous  year  [Figure  8(B)].  Clusters  classified  as  "isolated"  from 
other  active  clusters  under  the  k-means  classification  procedure  also  fledged 
fewer  young  than  clusters  that  were  relatively  less  isolated  [Figure  (8C)].  No 
consistent  differences  in  the  number  of  young  fledged  are  apparent  relative  to 
experience  of  the  breeding  male  or  female  in  the  cluster  [Figures  8(D)  and  8(E)], 


GRPN 


YEAR 

■  1997 
1998 

■  1999 


Figure  8  (A)  GRPN:  GRPN  “NO”  =  pairs  without  auxiliaries  and  “YES”  =  pairs  with  helpers. 
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PRSUC 

Figure  8  (B)  PRSUC:  PRSUC  “NO”  =  no  successful  nesting  attempts  in  the  previous  year  and  “YES”  = 
successful  nesting  attempt  in  the  previous  year. 


ISOL 


YEAR 

■  1997 
1998 

■  1999 


Figure  8  (C)  ISOL:  ISOL  “NO”  =  k-means  classification  categorical  indicating  reduced  cluster  isolation 
and  “YES”  =  k-means  classification  categorical  variable  indicating  increased  cluster  isolation. 


YEAR 

■  1997 
1998 

■  1999 


Figure  8  (D)  MBRYR:  See  Table  1  for  coding. 
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FBRYR 

Figure  8  (E)  FBRYR:  See  Table  1  for  coding. 


YEAR 

■  1997 
1998 

■  1999 


Figure  8.  Mean  number  of  fledged  by  year,  grouped  by  independent  variables. 


Table  6  summarizes  Mann-Whitney  U  test  statistics  and  probability  values  for 
means  of  number  of  young  fledged  grouped  by  independent  class  variables.  The 
dependent  variable,  mean  number  of  young  fledged,  was  significantly  different  in 
1998  and  1999  for  breeding  pairs  with  auxiliaries  versus  those  without  auxilia¬ 
ries.  Number  of  young  fledged  was  significantly  higher  in  all  years  in  clusters 
where  young  were  fledged  in  the  previous  year.  Although  clusters  that  were 
more  isolated  consistently  fledged  fewer  young  in  all  years,  the  means  were  only 
significantly  different  in  1999.  Neither  male  nor  female  breeding  tenure  in  the 
cluster  had  any  significant  effect  on  number  of  young  fledged  in  any  year;  how¬ 
ever,  sample  sizes  for  these  two  variables  were  relatively  small,  particularly  for 
adults  in  their  first  breeding  year. 


Table  6.  Non-parametric  means  test  for  NFLDG  grouped  by  demographic  variables. 


1997 

1998 

1999 

Grouping 

variables 

M-W  U* 

P 

N 

M-WU 

P 

N 

M-W  U 

P 

N 

GRPN 

78.5 

0.280 

29 

14.0 

0.000 

26 

77.5 

0.008 

36 

PRSUC 

48.5 

0.041 

31 

42.00 

0.041 

28 

74.5 

0.011 

37 

ISOL 

144.0 

0.530 

32 

130.5 

0.222 

29 

258.5 

0.045 

39 

MBRYR 

14.5 

0.736 

14 

25.5 

0.271 

23 

17.0 

0.215 

23 

FBRYR 

23.0 

0.452 

15 

20.0 

0.504 

15 

22.5 

0.380 

17 

*  M-W  U  =  Mann-Whitney  U  test  statistic,  p  =  significance  value  for  test  statistic.  N  =  total  sample  size. 
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The  other  dependent  variables  considered  in  these  analyses  were  class  variables 
with  two  levels.  Since  the  independent  variables  were  also  class  variables  with 
two  levels,  2x2  tables  were  analyzed  to  test  the  association  between  the  depend¬ 
ent  and  independent  variables.  Table  7  shows  Fisher's  exact  test  probabilities 
for  these  tests  of  independence.  Significance  values  are  not  shown  for  some  2x2 
tables  because  of  the  relationship  between  the  independent  and  dependent  vari¬ 
ables.  This  is  the  case  for  MBRYR  and  FBRYR  by  MALEREC  and  FEMREC, 
respectively.  In  this  case,  the  dependent  and  independent  variables  are  provid¬ 
ing  essentially  the  same  information  (breeding  tenure  in  the  cluster)  so  measure 
of  association  is  trivial.  A  similar  relationship  occurs  with  MALEREC  and 
FEMREC  with  PRSUC.  FEMREC  and  MALEREC  implies  in  most  cases  that  a 
breeding  male  was  not  present  in  the  previous  year,  which  also  implies  no  young 
fledged  in  the  previous  year.  The  variable  PRSUC  essentially  provides  the  same 
information.  There  is  also  a  trivial  comparison  for  GRPN  and  BREEDING, 
where  the  class  variable  GRPN  assumes  presence  of  a  breeding  pair.  I  n  several 
instances  there  were  no  observations  for  one  level  of  the  dependent  or  independ¬ 
ent  variable  so  the  2x2  table  could  not  be  constructed.  This  was  primarily  the 
case  for  those  variables  that  were  dependent  on  identification  of  banded  birds 
(dependent  variables  MALERET,  MALEREC,  FEMRET,  and  FEMREC,  and  in¬ 
dependent  variables  MBRYR  and  FBYR). 

Table  7  shows  that  the  only  consistently  significant  association  is  between 
whether  a  breeding  pair  is  present  and  whether  young  were  fledged  in  the  previ¬ 
ous  year  at  the  cluster.  No  other  tests  of  independence  were  consistently  re¬ 
jected  for  any  of  the  other  pairs  of  independent  and  dependent  variables.  Note, 
however,  that  sample  sizes  were  relatively  small  for  the  dependent  variables 
MALERET,  MALEREC,  FEMRET,  and  FEMREC,  and  the  independent  vari¬ 
ables  MBRYR  and  FBRYR. 

The  independent  variables  GRPN,  PRSUC  and  I  SOL  will  be  retained  for  evalua¬ 
tion  in  the  final  predictive  regression  models  for  NFLDG,  SUCCESS,  and 
BREEDING.  However,  the  dependent  variables  MALERET,  MALEREC, 
FEMRET,  and  FEMREC,  and  the  independent  variables  MBRYR  and  FBRYR 
will  not  be  evaluated  in  the  final  model.  These  variables  are  dependent  on  ob¬ 
servations  of  banded  birds  and  exhibit  a  high  frequency  of  missing  values  in  the 
data.  Further  analyses  of  these  variables  would  be  uninformative  due  to  low 
sample  sizes  resulting  from  missing  data. 
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Table  7.  Significance  values  for  Fisher’s  exact  test  of  independence  between  dependent 
variables  (rows)  and  independent  variables  (columns)  1997-1999. 


Independent  variables 

Dependent 

GRPN 

PRSUC 

ISOL 

MBRYR 

FBRYR 

variables 

BREEDING 

1997 

see  1 

0.000 

0.552 

0.250 

1.000 

1998 

0.000 

0.039 

1.000 

* 

1999 

0.000 

0.192 

* 

* 

SUCCESS 

1997 

0.354 

0.335 

1.000 

0.214 

1.000 

1998 

0.063 

0.142 

0.667 

1.000 

1.000 

1999 

0.104 

0.036 

0.273 

1.000 

* 

MALERET 

1997 

0.258 

1.000 

0.229 

see  2 

0.378 

1998 

1.000 

0.217 

1.000 

0.273 

1999 

0.093 

0.470 

0.560 

0.505 

MALEREC 

1997 

* 

see  3 

0.077 

* 

* 

1998 

* 

0.547 

* 

* 

1999 

* 

1.000 

* 

* 

FEMRET 

1997 

0.486 

0.008 

0.467 

0.378 

see  2 

1998 

0.228 

1.000 

1.000 

0.091 

1999 

0.620 

0.450 

1.000 

0.505 

FEMREC 

1997 

* 

see 3 

0.176 

* 

* 

1998 

* 

0.277 

* 

* 

1999 

* 

0.266 

* 

* 

*  No  values  observed  for  one  level  of  dependent  variable. 

1  GRPN  variable  implies  presence  of  breeding  pair.  Test  of  association  with  BREEDING  (presence  or 
absence  of  breeding  pair)  is  trivial. 

2  Dependent  and  independent  variables  are  essentially  providing  the  same  information;  e.g.,  a  male  in 
its  second  year  of  breeding  in  a  cluster  (MBRYR)  will  be  a  returning  breeder  from  the  previous  year 
(MALERET).  Tests  of  association  for  these  variables  are  trivial. 

3  MALEREC  and  FEMREC  implies  in  most  cases  that  a  breeding  male  was  not  present  in  the  previous 
year,  which  also  implies  no  young  fledged  in  the  previous  year.  The  variable  PRSUC  essentially  pro¬ 
vides  the  same  information  making  these  tests  of  independence  trivial. 


Habitat 

Before  analysis,  several  calculations  were  performed  on  collected  habitat  data. 
Basal  area  of  pines  and  hardwoods  within  clusters,  estimated  via  assessments 
with  a  ten-factor  prism,  was  converted  into  metric  units  (m2  per  hectare).  Pine 
basal  area  estimates  were  then  subdivided  into  two  size  classes  in  each  cluster: 
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trees  greater  than  25.4  cm  in  diameter  (overstory)  and  trees  less  than  25.4  cm  in 
diameter  (understory).  Number  of  stems  per  hectare  for  pines  and  hardwoods 
was  calculated  using  data  from  point-to-plant  measurements  taken  from  nest 
tree  centered  point-quarter  samples  in  each  cluster.  Stems  per  hectare  were  de¬ 
termined  for  each  point-quarter  plot,  and  the  average  of  these  measurements 
was  calculated  for  the  entire  cluster.  As  before,  pines  were  grouped  into  greater 
than  and  less  than  25.4-cm  size  classes.  Also,  because  variation  in  number  of 
stems  per  hectare  across  clusters  was  high,  estimates  were  log  transformed  be¬ 
fore  further  analysis. 

Percent  ground  cover  and  percent  canopy  cover  for  each  cluster  were  determined 
via  ocular  estimates  along  transects  in  each  of  the  five  point-quarter  plots.  On 
each  plot,  these  estimates  were  calculated  from  20  sightings  along  five  transects 
in  each  cluster.  Percentages  were  calculated  by  dividing  the  total  number  of 
positive  encounters  for  each  category  by  the  total  observations.  Both  canopy 
cover  and  ground  cover  were  subdivided  into  several  categories  (canopy  -  over¬ 
story  and  understory;  ground  -  woody  bushes,  grass  and  forbs,  open,  and  road). 
The  mean  of  the  estimates  from  each  of  the  five  point-quarter  plots  was  used  as 
an  index  of  the  cover  for  that  layer  for  the  entire  cluster. 

Because  many  of  these  habitat  variables  were  to  varying  degrees  correlated,  we 
reduced  the  total  number  of  variables  included  in  the  original  analysis  by  remov¬ 
ing  those  variables  that  showed  no  obvious  relationship  with  RCW  fitness  vari¬ 
ables.  Cover  data  were  not  collected  for  all  sample  clusters.  Analyses  of  avail¬ 
able  cover  data  indicated  no  relationships  between  cover  variables  and  RCW 
success  variables.  Therefore,  all  measures  of  cover  were  eliminated  from  further 
analysis  to  maintain  total  sample  size  for  other  variables  used  in  the  analysis. 

Of  the  variables  remaining  after  removal  of  cover  data,  several  were  strongly  col- 
linear  with  other  habitat  variables.  Because  inclusion  of  many  variables  in  a 
linear  regression  model,  regardless  of  the  strength  of  the  relationship  between 
predictive  and  predicted  variables,  artificially  elevates  the  predictive  power  of 
the  model,  it  was  necessary  to  remove  additional  variables.  Thus,  we  eliminated 
from  further  analysis  those  habitat  variables  exhibiting  strong  relationships  (r2 
>  0.70)  with  other  habitat  variables,  yet  relatively  weaker  relationships  with 
RCW  fitness  variables. 

Two  other  types  of  habitat-specific  data,  in  addition  to  the  previous  mentioned 
variables,  were  included  in  further  analysis:  cluster  burn  history  and  abundance 
of  cavities  per  cluster.  Burn  data  for  all  3  years  (1997-1999)  was  included,  but 
data  on  burn  history  are  somewhat  redundant.  Therefore,  only  data  from  1998, 
the  same  year  in  which  the  habitat  data  was  collected,  were  used  in  further 
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analysis.  Similarly,  only  data  from  1998  were  used  for  number  of  cavities  per 
cluster,  as  using  all  years  would  produce  problems  of  collinearity.  A  list  and  de¬ 
scription  of  all  habitat  variables  used  in  further  analysis  are  provided  in  Table  1 
(page  14). 

According  to  exploratory  correlation  analyses,  few  habitat  variables  exhibited 
obvious  relationships  with  RCW  fitness  variables.  We  therefore  used  principal 
component  analysis  (PCA),  a  multivariate  ordination  technique,  to  identify  the 
major  dimensions  of  variation  within  the  habitat  variable  matrix.  PCA  reduces 
dimensionality  of  a  data  matrix  by  combining  several  measured  variables  into  a 
few,  more  informative  variables,  somewhat  representative  of  a  "habitat  gradi¬ 
ent,"  and  identifies  those  measured  variables  that  are  most  influential  to  separa¬ 
tion  along  the  calculated  gradients  (see  Gauch  1982  for  a  detailed  explanation  of 
PCA).  All  listed  variables  were  included  in  this  ordination  procedure  except 
NCAV8.  This  variable  could  have  potentially  complicated  interpretation  of  gra¬ 
dients  suggested  by  PCA,  because  this  variable,  unlike  others,  is  not  directly  re¬ 
lated  to  the  habitat  itself,  but  is  a  combination  of  both  natural  conditions  and 
artificial,  anthropogenic  effects.  Though  burn  history  may  also  exhibit  an  an¬ 
thropogenic  aspect,  its  effects  on  habitat  should  be  quite  great,  whereas  cavity 
abundance  may  vary  irrespective  of  actual  habitat  gradients. 

To  perform  PCA,  we  first  relativized  all  included  variables  using  the  general  re- 
lativization  function  in  the  software  package,  PC-ORD  v3.20.  This  transforma¬ 
tion  maintains  variation  within  individual  variables,  yet  standardizes  variation 
between  variables,  therefore  eliminating  artifacts  in  ordination  due  to  non- 
uniform  variable  distributions  (McCune  1992).  Next,  we  used  the  vari¬ 
ance/covariance  PCA  procedure  to  ordinate  clusters  and  a  varimax  rotation  to 
maximize  separation  of  these  clusters  along  generated  principal  components 
(McCune  1992). 

Scores  from  the  first  three  principal  components  for  each  cluster  are  plotted  in 
three-dimensional  space  (Figure  9).  Table  8  provides  information  on  the  PCA 
axes  produced,  and  Table  9  provides  correlation  coefficients  between  PCA  scores 
and  raw  habitat  variables  after  varimax  rotation.  Axis  1,  which  explained  57.96 
percent  of  the  variation  in  the  main  habitat  matrix,  was  strongly  positively 
correlated  to  hardwood  basal  area  (Kendall's  tau  r  =0.891).  As  is  evident  in 
Figure  9,  most  clusters  along  this  axis  occurred  in  relatively  close  proximity  to 
each  other.  Very  few  clusters  had  high  or  low  hardwood  basal  areas  values;  most 
were  intermediate.  However,  though  it  appears  that  little  variation  exists  along 
this  axis,  this  variable  was  most  important  to  ordination  of  sites,  and  though 
only  a  few  clusters  occupied  extremes,  all  sites  were  tightly  arranged  along  this 
habitat  gradient  (see  above).  Thus,  based  on  the  evaluated  data  matrix, 
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between-site  variation  was  captured  by  this  axis.  Axis  2  explained  20.84  percent 
of  the  overall  variation  in  the  habitat  matrix  and  was  strongly  positively 
correlated  to  understory  pine  basal  area  (Kendall's  tau  r  =0.910).  This  axis, 
though  explaining  less  variation  than  the  first,  better  separated  RCW  clusters 
along  its  length.  Midstory  pine  basal  area  was  more  evenly  distributed  among 
sites  than  was  hardwood  basal  area.  Axis  3  explained  12.52  percent  of  the 
variation  in  the  main  matrix  and  was  strongly  negatively  correlated  to  burn 
interval  1998  and  positively  correlated  to  number  of  burns  in  the  past  5  years 
(Kendall's  tau  r  =  -  0.766,  and  0.573,  respectively).  RCW  clusters  were  well 
spread  evenly  along  this  component.  An  additional  4  axes  explained  the 
remaining  8.7  percent  variation  in  the  main  matrix;  however,  because  the 
explanatory  power  of  these  components  was  small  in  comparison  to  the  first 
three  axes  and  because  no  clear  relationship  between  PC  scores  and  measured 
habitat  variables  were  identified  with  these  axes,  components  four  through  six 
were  eliminated  from  further  consideration. 


PC2 


PCI  positively  correlated  to  increasing  hardwood  basal  area,  PC2  positively  correlated  with  increasing 
midstory  pine  basal  area,  and  PC3  negatively  correlated  to  1998  burn  interval. 


Figure  9.  Ordination  of  RCW  clusters  along  first  three  principal  components  of  measured 
habitat  variables. 
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Table  8.  Critical  values  from  principal  component  analysis  of  RCW  cluster  habitat  variables. 


Axis 

Eigenvalue 

%  of  Variance 

Cumulative  %  of  Variance 

1 

0.033 

57.955 

57.955 

2 

0.012 

20.841 

78.796 

3 

0.007 

12.516 

91.312 

4 

0.002 

3.657 

94.968 

5 

0.002 

2.430 

97.398 

6 

0.001 

1.841 

99.240 

7 

0.001 

0.760 

100.000 

Table  9.  Correlation  coefficients  of  PCA  scores  and  measured  RCW  cluster  habitat  variables 
after  varimax  rotation. 


Kendall’s  Tau  correlation  coefficient  with  PCA  axis 

Habitat  variable 

1 

2 

3 

PINEBA 

0.080 

0.280 

0.000 

PINEL10 

-0.163 

0.910 

-0.161 

HDBA 

0.891 

-0.002 

-0.207 

PINESTML 

0.007 

0.145 

-0.133 

HDSTML 

-0.026 

0.296 

0.225 

5BURN8 

-0.491 

0.021 

0.573 

IBURN8 

0.407 

0.037 

-0.766 

To  assess  univariate  relationships  between  habitat  characteristics  and  RCW 
success  measures,  two  statistical  procedures  were  used.  First,  Spearman  rank 
correlations  between  PC  scores  and  number  of  cavities  in  1998  and  number  of 
RCW  fledglings  per  cluster  (NFLDG)  in  all  3  years  were  calculated  (Table  10). 
No  significant  relationships  were  detected  in  these  comparisons.  Furthermore, 
not  only  was  significance  absent,  but  also  relationships  between  habitat  vari¬ 
ables  and  number  of  fledglings  across  years  were  in  opposition.  For  example, 
PC2  was  positively  related  to  RCW  success  in  1997  but  negatively  associated  in 
1998  and  1999.  This  inconsistency  suggests  that  the  lack  of  a  relationship  be¬ 
tween  habitat  variables  and  number  of  fledglings  cannot  be  explained  solely  by 
low  sample  size. 
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Table  10.  Spearman  rank  correlation  coefficients  between  principal  components  scores  of 
measured  habitat  variables  and  number  of  useable  RCW  cavities  in  1998  and  number  of 
fledglings  produced  by  cluster  for  years  1997, 1998,  and  1999. 


Habitat  Variable* 

#  Fledglings  1997 

#  Fledglings  1998 

#  Fledglings  1999 

PCI 

-0.349 

-0.256 

-0.104 

PC2 

0.046 

-0.037 

-0.058 

PC3 

0.127 

-0.238 

-0.225 

#  cavities  1 998 

-0.027 

0.328 

0.198 

*  No  relationships  were  statistically  significant  (P  <  0.05). 


Next,  logistic  regression  was  used  to  assess  relationships  between  habitat  vari¬ 
ables  and  binomial  RCW  success  variables  including  breeding  pair  presence 
(BREEDING),  successful  fledging  (SUCCESS),  male/female  recruitment 
(MREC/FREC),  and  male/female  retention  (MRET/FRET)  (Table  11).  As  with 
fledgling  number,  no  consistent  relationships  were  identified  across  all  three 
years,  though  significant  predictive  power  was  indicated  in  two  cases:  breeding 
pair  presence  in  1997  by  PC2  scores  and  in  1999  by  number  of  cavities.  Such  re¬ 
lationships,  however,  are  likely  spurious. 

In  conclusion,  it  appears  that  though  principal  component  analysis  adequately 
captured  habitat  variability,  this  variability  did  not  explain  variation  in  RCW 
success  across  clusters:  no  significant  trends  between  principal  component  scores 
and  measures  or  RCW  success  were  detected.  F  urthermore,  the  number  of  cavi¬ 
ties  present  in  the  area,  a  variable  that  has  been  shown  previously  to  be  influen¬ 
tial  to  RCW  success,  was  not  significantly  predictive  of  RCW  fitness  here. 

Such  lack  of  relationships  may  be  due  to  several  factors.  First,  it  is  possiblethat 
in  fact  no  relationships  do  exist.  Most  previous  studies  of  RCW  success  have  had 
difficulty  identifying  specific  habitat  variables  that  are  predictive  of  RCW  suc¬ 
cess.  Most  variables  that  have  been  identified  are  not  large  scale,  physiognomic 
variables  like  the  ones  identified  here,  but  instead  they  are  nest  site  specific. 
Second,  habitat  conditions  among  clusters  evaluated  in  this  study  may  not  have 
been  sufficiently  different  to  affect  RCW  success.  Though  PCA  did  separate  clus¬ 
ters  in  multidimensional  space,  it  is  possiblethat  little  variation  in  quality  as 
assessed  by  RCWs  occurred  along  these  axes.  Lastly,  variables  measured  and 
included  in  the  analysis  may  not  have  been  those  used  by  RCWs  as  measures  of 
"good"  habitat.  As  apparent  from  the  data,  certain  clusters  were  consistently 
more  successful  across  years  than  others.  This  trend  implies  that  certain  char¬ 
acteristics  of  the  area  do  in  fact  influence  RCW  success.  Perhaps  the  variables 
measured  did  not  identify  these  characteristics,  thus  no  significant  habitat-RCW 
success  relationships  were  indicated. 
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Table  11.  Critical  values  (Wald  statistic  (W)  and  P-value  (P))  from  logistic  regression  analyses  of 
principal  component  scores  of  habitat  variables  and  number  of  useable  RCW  cavities  in  1998 
and  binomial  measures  of  RCW  reproductive  success  in  the  years  1997-1999. 


Success  measure  Habitat  variables 

PCI  PC2  PC3  #  cavities98 


W 

P 

W 

P 

W 

P 

W 

P 

Nest  success 

1997 

3.46 

0.06 

0.67 

0.41 

1.77 

0.18 

1.06 

0.30 

1998 

1.38 

0.24 

1.88 

0.17 

0.01 

0.92 

2.68 

0.10 

1999 

1.02 

0.31 

0.22 

0.64 

0.53 

0.47 

0.03 

0.86 

Breeding  pair 

presence 

1997 

0.02 

0.89 

6.00 

0.01* 

0.02 

0.90 

2.43 

0.12 

1998 

1.09 

0.30 

0.64 

0.42 

0.41 

0.52 

2.46 

0.12 

1999 

0.18 

0.67 

0.21 

0.65 

2.90 

0.88 

5.01 

0.03* 

Male  retention 

1997 

1.99 

0.16 

0.16 

0.69 

0.00 

0.97 

1.86 

0.17 

1998 

sample  size  insufficient 

1999 

0.91 

0.34 

1.71 

0.19 

0.27 

0.61 

0.42 

0.52 

Female  retention 

1997 

sample  size  insufficient 

1998 

0.01 

0.92 

0.64 

0.42 

0.67 

0.41 

0.04 

0.83 

1999 

0.34 

0.56 

0.36 

0.55 

0.79 

0.38 

0.16 

0.69 

Male  recruitment 

1997 

1.80 

0.18 

0.97 

0.33 

0.38 

0.54 

2.55 

0.11 

1998 

0.68 

0.41 

2.14 

0.14 

0.21 

0.65 

1.33 

0.25 

1999 

1.29 

0.26 

2.19 

0.14 

0.70 

0.40 

0.13 

0.72 

Female  recruitment 

1997 

0.11 

0.74 

0.00 

0.95 

0.09 

0.77 

0.04 

0.83 

1998 

0.00 

0.97 

1.01 

0.31 

0.33 

0.56 

1.13 

0.29 

1999 

0.03 

0.86 

0.02 

0.90 

0.57 

0.45 

0.57 

0.45 

indicates  significant  P-value 

(<  0.05) 

Human  Disturbance 

Measures  of  human  disturbance  in  proximity  to  RCW  clusters  were  examined  for 
effects  on  RCW  demographic  measures.  Dependent  variables  in  these  analyses 
were  nest  success  in  clusters  with  a  breeding  pair  (SUCCESS;  did  a  nest  fledge 
at  least  one  young),  number  of  young  fledged  from  clusters  with  a  breeding  pair 
(NFLDG),  presence  of  breeding  pair  (BREEDING),  male  and  female  adult 
retention  (MALERET,  FEMRET;  retention  of  the  breeding  adult  from  the 
preceding  year)  and  male  and  female  adult  recruitment  (MALEREC,  FEMREC). 
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Independent  (predictor)  variables  used  were  the  proportion  of  10-minute 
observations  in  which  military  activity  (presence  or  transit  of  military  vehicles  or 
personnel)  was  observed  (MILACT,  arcsine  transformation),  proportion  of  10- 
minute  observations  in  which  any  human  activity  (military  and  civilian  activity) 
was  observed  (ACT,  arcsine  transformed),  distance  to  maintained  roads 
(LIOROAD,  log  transformed)  and  number  of  scheduled  training  days  (N  DAY). 

Table  12  shows  Spearman's  rank  correlation  matrices  for  number  of  young 
fledged  (NFLDG),  total  activity  observed  (ACT),  military  activity  observed 
(MILACT)  and  distance  to  maintained  roads  for  the  years  1997-1999  (ROAD). 
No  strong  correlations  were  found  between  number  of  young  fledged  and  any  of 
the  independent  measures  of  potential  disturbance.  Total  activity  was  highly 
correlated  (Spearman’s  rank  correlation  coefficient  >0.750)  with  military  activ¬ 
ity  in  all  years.  Due  to  the  strength  of  this  correlation,  the  variable  and  deriva¬ 
tives  of  observed  military  activity  were  used  for  all  subsequent  analyses. 


Table  12.  Spearman  rank  correlation  matrix  for  potential  disturbance  variables  and  number  of 
young  fledged,  1997-99. 


Year 

NFLDG 

ACT* 

MILACT 

ROAD 

1997 

NFLDG 

1.000 

ACT 

-0.430 

1.000 

MILACT 

-0.262 

0.837 

1.000 

ROAD 

0.162 

-0.590 

-0.832 

1.000 

1998 

NFLDG 

1.000 

ACT 

0.125 

1.000 

MILACT 

0.150 

0.805 

1.000 

ROAD 

-0.009 

-0.281 

-0.367 

1.000 

1999 

NFLDG 

1.000 

ACT 

0.123 

1.000 

MILACT 

0.060 

0.759 

1.000 

ROAD 

0.066 

-0.288 

-0.335 

1.000 

*  ACT  and  MILACT  are  proportions  of  all  activity  and  military  activity  observed  during  10-minute  sample 
periods.  Data  are  arcsine  transformed.  ROAD  is  log  transformed  distance  of  nest  cavity  tree  or  nearest 
active  cavity  tree  to  a  maintained  road.  NFLDG  is  number  of  young  fledged  per  breeding  pair. 


Table  13  shows  results  of  non-parametric  Kruskal-Wallis  two-sample  test  mean 
number  fledged  on  factors  SITE  LOAD  and  ROAD(CLASS).  SITE  LOAD  and 
ROAD(CLASS)  are  class  variables  with  two  levels  derived  by  k-means 
classification  from  MILACT  and  ROAD.  These  univariate  means  tests  showed 
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no  statistically  significant  difference  (p  >0.05)  in  any  year  in  number  of  young 
fledged  from  clusters  with  higher  levels  of  observed  activity  versus  those  with 
low  levels  of  observed  military  activity  and  no  difference  in  number  of  young 
fledged  from  clusters  close  to  roads  versus  those  distant  from  roads.  Figures  10 
and  11  show  mean  number  of  young  fledged  by  year  for  these  two  factors.  Mean 
number  fledged  was  consistently  slightly  lower  for  clusters  closer  to  roads  for  all 
years  (Figure  11)  but  only  slightly  so.  No  consistent  trend  was  observed  for  the 
factor  SITE  LOAD. 


Table  13.  Non-parametric  means  test  for  NFLDG  grouped  by  SITELOAD  and  ROAD  (CLASS). 


SITELOAD* 

ROAD  (CLASS) 

Mann- 
Whitney  U 

P 

N 

Mann- 
Whitney  U 

P 

N 

NFLDG 

1997 

31.0 

0.469 

19 

127.0 

0.472 

32 

1998 

98.0 

0.473 

29 

78.0 

0.610 

29 

1999 

133.5 

0.731 

39 

181.0 

0.229 

39 

*  SITELOAD  and  ROAD  (CLASS)  variables  have  two  levels  based  on  k-means  classification.  N  = 
sample  size.  P  =  significance  values  for  Mann-Whitney  U  test  statistic. 


SITELOAD 


YEAR 

■  1997 
1998 

■  1999 


Figure  10.  Mean  number  of  young  fledged  by  year,  grouped  by  SITELOAD 
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YEAR 

■  1997 

■  1998 

■  1999 

FAR  NEAR 
ROAD(CLASS) 

Figure  11.  Mean  number  of  young  fledged  by  year,  grouped  by  ROAD(CLASS). 

The  other  dependent  variables  considered  in  these  analyses  are  class  variables 
with  two-levels.  Logistic  regressions  were  performed  with  MILACT  and  ROAD 
as  predictor  variables  for  each  of  the  independent  variables  SUCCESS, 
BREEDING,  MALERET,  MALEREC,  FEMRET,  and  FEMREC.  Results  of  these 
logistic  regressions  for  all  years  are  presented  in  Table  14.  The  independent 
variables  were  poor  predictors  for  all  the  dependent  variables  in  all  years  (p- 
values  for  logit  coefficients  >0.05).  Overall  model  fit  was  also  generally  poor  as 
indicated  by  the  lack  of  significance  for  the  log-likelihood  ratio.  The  only  signifi¬ 
cant  log-likelihood  ratio  was  for  SUCCESS  in  1998  (p  =0.047)  and  FEMRET  in 
1999  (p  =0.003);  however,  examination  of  prediction  success  tables  for  these  two 
regressions  indicated  the  model  was  a  poor  predictor  of  these  dependent  vari¬ 
ables.  Also,  model  fit  for  these  dependent  variables  in  alternate  years  did  not 
corroborate  the  model  fit  for  these  dependent  variables.  Note  that  sample  sizes 
are  low  in  many  cases  for  MALERET,  MALEREC,  FEMRET,  and  FEMREC,  and 
consequently  the  regression  statistics  are  suspect. 
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Table  14.  Logistic  regression  statistics. 


Independent  Variables 

MILACT 

J3 

o 

> 

D 

CLASS) 

Dependent 

Variables* 

t-ratio 

p-value 

t-ratio 

p-value 

Log-likelihood 
ratio  (p-value) 

N 

1997 

SUCCESS 

0.015 

0.988 

0.049 

0.961 

0.004  (0.998) 

19 

BREEDING 

-0.587 

0.557 

-1.847 

0.065 

5.467  (0.065) 

24 

MALERET 

0.031 

0.976 

0.415 

0.678 

1.976  (0.372)** 

9 

MALEREC 

-0.006 

0.995 

-0.586 

0.558 

1.421  (0.492)* 

4 

FEMRET 

0.030 

0.976 

0.457 

0.648 

1.100  (0.577)** 

12 

FEMREC 

0.016 

0.988 

-0.023 

0.982 

4.499  (0.105)** 

4 

1998 

SUCCESS 

1.621 

0.105 

-0.380 

0.704 

6.120(0.047) 

29 

BREEDING 

0.953 

0.341 

1.784 

0.074 

3.965  (0.138) 

51 

MALERET 

0.032 

0.974 

-0.352 

0.725 

1.560  (0.458)** 

24 

MALEREC 

-0.986 

0.324 

0.933 

0.351 

2.330  (0.312) 

12 

FEMRET 

0.034 

0.973 

-0.002 

0.999 

7.191  (0.027)** 

14 

FEMREC 

-1.343 

0.179 

0.920 

0.358 

3.364  (0.186) 

14 

1999 

SUCCESS 

-0.027 

0.978 

0.851 

0.395 

0.794  (0.627) 

39 

BREEDING 

0.950 

0.342 

1.302 

0.193 

2.563  (0.278) 

51 

MALERET 

-1.193 

0.233 

-0.333 

0.739 

1.564  (0.458) 

23 

MALEREC 

1.461 

0.140 

0.525 

0.600 

4.067  (0.131) 

12 

FEMRET 

0.582 

0.560 

-1.740 

0.082 

11.337  (0.003) 

17 

FEMREC 

-0.544 

0.586 

1.181 

0.238 

3.474  (0.176) 

14 

*  Dependent  variables  =  Constant  +  MILACT  +  ROAD.  SUCCESS  =  fledging  at  least  one  young  given 
presence  of  breeding  pair.  BREEDING  =  presence  of  breeding  pair.  MALERET  =  presence  of  breeding 
male  from  previous  year.  MALEREC  =  recruitment  of  new  breeding  male.  FEMRET  and  FEMREC  = 
female  equivalent  of  MALERET  and  MALEREC  variables.  MILACT  and  ROAD  are  arcsine  and  log- 
transformed  values,  respectively.  N  =  sample  size. 

**Model  failed  to  satisfy  change  tolerance.  Log-likelihood  ratio  and  associated  p-value  are  suspect. 
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6  Analyses  of  Maneuver  Training  Effects 

Chapter  5  presented  analyses  and  data  summaries  to  select  variables  for  analy¬ 
sis  of  maneuver  training  effects  on  RCWs  relative  to  habitat  and  demographic 
factors.  Variables  were  retained  based  on  correlations  among  independent  vari¬ 
ables,  relationship  to  dependent  variables,  sample  size,  and  representative  of 
factors  potentially  affecting  RCW  populations.  Dependent  variables  retained  for 
analysis  include  number  of  young  fledged  (NFLDG,  Table  1  (page  14),  whether 
young  were  fledged  from  the  cluster  (SUCCESS),  and  whether  a  breeding  pair 
was  present  (BREEDING).  Independent  variables  retained  were  presence  of 
auxiliary  birds  in  breeding  groups  (GRPN),  success  in  fledging  young  in  the  pre¬ 
ceding  year  (PRSUC),  principal  component  scores  for  habitat  (PCI,  PC2,  and 
PC3),  level  of  military  activity  observed  at  cluster  site  (Ml  LACT),  local  density  of 
active  clusters  (I  SOL),  number  of  useable  cavities  per  cluster  (NCAV),  and  prox¬ 
imity  of  cluster  to  nearest  maintained  road  (LIOROAD). 

Analyses  of  the  relationships  between  dependent  and  independent  variables 
were  performed  using  multiple  linear  regression  approaches  (GLM  and  logistic). 
In  a  recent  review,  MacNally  (2000)  discusses  the  potential  pitfalls  in  ecological 
applications  of  this  approach,  particularly  in  predictive  versus  causal  inferences 
of  the  results.  The  results  reported  here  present  the  most  parsimonious  predic¬ 
tive  relationships  between  the  dependent  and  independent  variables.  These  re¬ 
lationships  are  used  to  evaluate  the  relative  importance  of  the  independent  vari¬ 
ables  as  predictors  of  observed  variance  in  the  independent  variables. 

The  dependent  variables  SUCCESS  and  BREEDING  are  class  variables  with 
two  levels  (Yes  or  No).  Logistic  regression  was  used  to  evaluate  predictive  mod¬ 
els  for  SUCCESS  and  BREEDING.  Results  of  model  exploration  were  satisfac¬ 
tory  in  identifying  any  combination  of  significant  predictor  variables  for 
SUCCESS  or  BREEDI NG.  For  both  dependent  variables,  models  either  failed  to 
satisfy  change  tolerance,  or  in  stepwise  procedures  (backward  elimination)  all 
independent  variables  except  the  constant  were  removed  due  to  lack  of  signifi¬ 
cance.  These  results  suggest  that  the  independent  variables  or  combination  of 
variables  are  not  reliable  predictors  of  the  presence  of  a  breeding  pair  or  success 
of  breeding  pairs  in  fledging  young. 

Multiple  regression  to  evaluate  predictive  models  for  number  of  young  fledged 
(NFLDG)  produced  more  satisfying  results.  Table  15  shows  results  for  multiple 
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linear  regression  for  the  full  model  (all  independent  variables  entered).  R2  for 
the  full  model  ranged  from  0.404  in  1999  to  0.723  in  1998.  The  regression  sum- 
of-squares  was  significant  for  1998  (F  =4.361,  p  =0.006).  However,  as  MacNally 
(2000)  points  out,  with  many  independent  variables  in  the  model,  high  R2  values 
can  be  a  mathematical  artifact  of  the  proportion  of  the  number  of  variables  to 
cases.  Also  the  purpose  is  to  identify  those  variables  that  are  the  best  predictors 
of  the  number  of  young  fledged. 


Table  15.  Results  of  multiple  linear  regression  for  NFLDG  for  all  independent  variables  entered 
as  model  terms. 


Year 

R2 

F 

P 

1997 

0.503 

0.674 

0.715 

1998 

0.723 

4.361 

0.006 

1999 

0.404 

1.816 

0.117 

To  evaluate  the  most  parsimonious  model  for  prediction  of  number  of  young 
fledged,  a  stepwise  regression  (backward  elimination)  was  performed.  Criteria 
for  variable  removal  was  p  =0.15.  Again  MacNally  (2000)  cautions  about  the 
potential  to  retain  "spurious”  variables  based  on  p-values.  However,  in  this 
analysis,  model  results  from  different  years  can  be  considered  an  independent 
cross-validation  of  model  results. 

Table  16  shows  results  of  stepwise  regression  for  each  of  the  years  in  this  study. 
For  each  year  variables  are  ordered,  based  on  order  of  removal  and  significance 
values  of  final  variables  retained.  Note  that  although  the  stepwise  procedure 
converges  on  slightly  different  combinations  of  predictor  variables  depending  on 
the  year,  the  relative  order  in  which  variables  are  retained  is  fairly  consistent 
among  years,  particularly  among  those  most  likely  to  be  retained  in  the  final 
variables.  GRPN,  PC3,  and  PRSUC  were  among  the  top  three  variables  in  at 
least  2  of  the  3  years.  The  only  other  variable  occurring  among  the  top  three  was 
I  SOL  in  one  year,  1998.  PRSUC  in  1998  was  fourth  .  Also  note  that  MILACT 
was  in  the  bottom  three  (variables  with  least  significance)  in  all  years. 
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Table  16.  Comparison  of  relative  significance  of  regression  coefficients  of  independent 
variables  in  stepwise  regression  (GLM)  for  number  of  young  fledged. 


1997* 

1998 

1999 

Most  significant 

PRSUC** 

GRPN** 

GRPN** 

GRPN 

PC3** 

PC3** 

PC3 

ISOL** 

PRSUC** 

LIOROAD 

PRSUC 

PC2 

PC2 

PC2 

NCAV 

ISOL 

PCI 

ISOL 

MILACT 

NCAV 

MILACT 

PCI 

MILACT 

PCI 

Least  significant 

NCAV 

LIOROAD 

LIOROAD 

*Relative  significance  is  shown  for  all  3  years.  For  each  year,  variables  are  ordered  from  most  signifi¬ 
cant  regression  coefficient  (more  likely  to  be  retained  in  stepwise  procedure)  to  least  significant  (most 
likely  to  be  removed  in  stepwise  procedure). 

**  Identifies  variables  retained  after  completion  of  stepwise  elimination.  Criterion  for  variable  elimination 
was  p  =  0.15. 


The  most  parsimonious  model  will  identify  the  fewest  model  terms  that  will  sig¬ 
nificantly  and  most  consistently  account  for  observed  variance  in  the  number  of 
young  fledged.  Based  on  the  results  above,  models  with  the  independent  vari¬ 
ables  GRPN,  PRSUC,  and  PC3  were  compared  for  the  full  model  including  all 
model  terms  and  all  subsets  of  the  full  model  with  two  model  terms.  Table  17 
shows  results  for  these  models.  Regression  sum-of-squares  was  significant  in  all 
years  for  the  models  GRPN  +  PRSUC  and  PC3  +  PRSUC.  The  regression  sum- 
of-squares  was  significant  in  1998  and  1999  for  the  full  model  and  GRPN  +PC3, 
but  were  not  significant  in  1997,  although  the  full  model  approached  significance 
(p  =0.073). 


Table  17.  Comparison  between  complete  multivariate  regression  model  (GLM)*  results  for 
NFLDG  with  all  possible  submodels  with  two  independent  variables. 


1997 

1998 

1999 

Model 

Terms 

R2 

F 

P 

R2 

F 

P 

R2 

F 

P 

GRPN+ 

PC3+ 

PRSUC 

0.247 

2.631 

0.073 

0.554 

8.700 

0.001 

0.364 

5.723 

0.003 

GRPN+ 

PC3 

0.022 

0.299 

0.744 

0.566 

14.977 

0.000 

0.303 

7.156 

0.003 

GRPN+ 

PRSUC 

0.219 

3.503 

0.046 

0.439 

8.615 

0.002 

0.294 

6.470 

0.004 

PC3+ 

PRSUC 

0.298 

5.953 

0.007 

0.247 

4.096 

0.029 

0.261 

5.989 

0.006 

*  Complete  GLM  model  is  NFLDG  =  constant  +  GRPN  +  PC3  +  PRSUC. 
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Examination  of  R2  values  for  the  full  model  and  submodels,  shows  that  the  full 
model  had  consistently  higher  R2  values  in  all  years  except  PC3  +  PRSUC  in 

1997  (R2  =  0.298  versus  0.247  for  full  model)  and  slightly  higher  values  for 
GRPN+PC3  in  1998  (R2  =  0.566  versus  0.554  for  full  model).  The  results  pre¬ 
sented  in  Table  17  suggest  that  of  all  the  predictor  variables  considered  in  this 
chapter  the  model  NFLDG  =GRPN  +  PC3  +  PRSUC  presents  the  best  balance 
between  accounting  for  the  largest  amount  of  variance  in  NFLDG  and  achieving 
significance  for  the  regression.  Significance  of  the  model  terms  GRPN  and 
PRSUC  for  NFLDG  is  consistent  with  univariate  results  presented  in  Chapter  5. 
It  is  also  of  interest  that  PC3  loads  most  heavily  on  variables  related  to  burning 
regime.  Specifically,  the  results  here  suggest  a  relationship  of  NFLDG  with 
higher  burn  frequency  and  interval  to  most  recent  burn  (as  suggested  by  PC3). 

Note,  however,  that  this  three-term  model  does  not  account  for  a  large  propor¬ 
tion  of  the  variance  in  NFLDG.  The  model  terms  retained  in  this  model  (GRPN, 
PC3,  and  PRSUC)  also  reinforces  that  these  models  are  predictive  and  do  not  il¬ 
lustrate  causal  relationships.  It  is  intuitive  that  if  a  cluster  is  successful  in 
fledging  young  in  the  previous  year  (PRSUC)  that  the  likelihood  of  having  auxil¬ 
iaries  (GRPN)  increases  in  the  following  year.  The  fact  that  young  were  fledged 
in  a  previous  year  should  not  be  interpreted  as  the  cause  of  increased  production 
of  young  in  the  current  year.  Rather  the  predictive  relationship  of  PRSUC  to 
NFLDG  likely  reflects  a  relationship  of  PRSUC  to  the  presence  of  auxiliaries  in 
the  current  year  (GRPN)  and  longer  tenure  of  the  breeding  adults  in  the  cluster. 
Although  the  relationship  is  not  exceedingly  strong,  correlation  between  PRSUC 
and  GRPN  are  among  the  highest  observed  among  the  independent  variables  in 

1998  (Spearman  rank  coefficient  =0.446)  and  1999  (Spearman  rank  coefficient  = 
0.502). 
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7  Determining  Population  Level  Effects  of 
Military  Disturbance 

Although  only  a  few  RCW  clusters  on  Fort  Stewart  appear  to  be  subject  to  rela¬ 
tively  high  levels  of  human  activity  at  present,  this  situation  could  change  after 
the  implementation  of  the  revised  1996  "Management  Guidelines  for  RCWs  on 
Army  Installations,”  which  are  anticipated  to  be  in  effect  by  2001.  The  question 
that  will  be  addressed  in  this  section  is:  What  hypothetical  effects  on  Red- 
cockaded  Woodpecker  (RCW)  population  viability  would  accrue,  if  various  pro¬ 
portions  of  RCW  breeding  habitat  on  Fort  Stewart  were  subject  to  high  levels  of 
potentially  disturbing  human  activity?  An  answer  to  this  question  will  be  at¬ 
tempted  by  performing  a  population  viability  analysis,  using  the  data  on  mili¬ 
tary-civilian  activities,  and  on  RCW  nest  productivity,  that  were  collected  at  Fort 
Stewart  RCW  cluster  sites  during  this  study. 


Materials  and  Methods 

The  procedure  for  estimating  the  potential  effects  of  military-civilian  activities 
on  the  Fort  Stewart  RCW  population  was  as  discussed  in  the  following  para¬ 
graphs. 

Classification  of  Cluster  Sites  by  Military  Activity  Levels 

First,  51  cluster  sites  having  at  least  2  years  of  observation  data  on  military  and 
civilian  activity  were  identified.  From  these,  a  subset  of  24  clusters,  each  with  3 
years  of  breeding  attempts  (breeding  pair  with  nest  present)  recorded  per  clus¬ 
ter,  was  established  as  the  data  set  for  the  present  analysis.  This  data  set  was 
expressed  as  a  bivariate  plot,  with  each  point  representing  a  cluster's  2-  or  3- 
year  (depending  on  available  data)  mean  proportion  of  observation  periods  with 
civilian  (mean  CIVACT)  vs.  military  (mean  MILACT)  activity.  For  definitions  of 
MILACT  and  CIVACT,  see  Table  1  (page  14).  Thedusters  were  then  divided,  on 
the  basis  of  the  bivariate  plot,  into  classes  of  "high”  vs.  "low”  military-civilian  ac¬ 
tivity  potentially  disturbing  to  RCWs  on  Fort  Stewart.  "High"  activity  classes 
included  those  clusters  where  the  proportion  of  observation  periods  with  civilian 
activity  was  greater  than  0.10,  and  the  proportion  of  observation  periods  with 
observed  military  activity  was  greater  than  0.15. 
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Fledgling  Production  in  Each  Activity  Level 

The  yearly  descriptive  statistics  for  "nest  success”  (proportion  of  nests  for  which 
SUCCESS  =1),  and  for  the  mean  number  of  fledglings  produced  (mean  NFLGD) 
per  successful  nest  (denoted  here  as  "Fledglings  |  Successful"),  were  calculated 
for  the  "low"  and  "high"  disturbance  activity  classes.  For  definitions  of  NFLGD 
and  SUCCESS,  see  Table  1.  The  means  of  these  yearly  statistics  were  tested  for 
difference  between  activity  classes  by  t-tests,  using  SAS/STATa  PROC  TTEST 
(SAS  I  nstitute  I  nc.  1989). 

Next,  the  means  for  nest  success  and  fledglings  per  successful  nest,  along  with 
their  temporal  (over  3  years)  and  sampling  error  variances  and  standard  devia¬ 
tions,  were  calculated  for  each  disturbance  activity  class  of  RCW  clusters.  Tem¬ 
poral  and  sampling  variances  (and  standard  deviations)  of  these  variables  were 
calculated  using  the  variance  partitioning  method  described  in  Stewart-Oaten, 
Murdoch,  and  Walde  (1995).  These  statistics  were  taken  to  represent  the  popu¬ 
lation  in  pure  "low"  vs.  pure  "high"  disturbance  activity  habitats.  Effective 
sample  sizes  for  the  low  and  high  activity  groups  were  taken  to  be  the  total 
number  of  observations  upon  which  the  analyses  were  based.  These  were  21 
clusters  x  3  years  =  63  cluster-years  for  the  "low"  activity  class,  and  a  corre¬ 
sponding  3  clusters  x  3  years  =9  cluster -years  for  the  "high"  activity  class. 

Hypothetical  Disturbance  Scenarios 

Weighted  averages  of  the  statistics  from  these  "pure"  activity  classes  were  used 
to  estimate  nest  productivity  statistics  (and  their  effective  sample  sizes)  for  nine 
hypothetical  disturbance  scenarios,  representing  reproductive  productivity  on 
Fort  Stewart  if  different  proportions  of  the  RCW  habitat  were  in  the  "high"  dis¬ 
turbance  activity  category.  These  scenarios  allow  comparison  of  the  observed 
disturbance  levels  on  Fort  Stewart  in  the  context  of  the  full  range  of  possibledis- 
turbance  scenarios.  For  given  statistics  from  "high"  vs.  "low"  activity  classes, 
Xhigh  and  Xiow,  the  weighted  average  for  the  scenario  with  proportion  P  of  the 
habitat  in  the  "high"  activity  category  was  calculated  as: 

Weighted  Average  of  X  =P  Xhigh  +  (1-P)  Xiow 

Effective  sample  sizes  for  each  scenario  were  also  calculated  as  weighted  aver¬ 
ages,  truncated  to  their  integer  values. 
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Population  Viability  Analysis  Using  PVAvES  1.0 

Finally,  the  statistics  for  these  nine  scenarios  were  entered  as  input  data  to  pro¬ 
gram  PVAvES  1.0,  a  population  viability  program  developed  at  the  U.S.  Army 
Engineer  Research  and  Development  Center /Construction  Engineering  Research 
Laboratory  (ERDC/CERL)  to  estimate  the  risk  of  extinction  for  threatened  and 
endangered  avian  populations  on  military  lands.  The  algorithms  of  the  PVAvES 
program  are  described  in  a  previous  technical  report  (Melton  et  al.  2001).  Deri¬ 
vations  of  many  of  the  input  parameter  values  used  in  the  present  study,  calcu¬ 
lated  from  Fort  Stewart  data,  are  described  there.  I  nput  parameter  derivations 
from  data,  other  than  those  noted  above,  will  not  be  described  in  detail  unless 
the  values  used  here  differed  from  those  of  the  previous  report. 

The  Appendix  shows  an  example  of  an  input  data  file  read  by  PVAvES  at  the  be¬ 
ginning  of  program  operation.  Parameter  values  from  the  Fort  Stewart  RCW 
population  used  in  this  study  are  as  follows: 

1.  Parameters  with  the  same  values  in  all  disturbance  scenarios: 

a.  The  initial  population  size  (Appendix,  input  #7)  is  set  here  to  the  known 
breeding  female  population  size  at  the  end  of  year  1999,  which  was  165 
breeding  females  (i.e.,  females  known  to  have  made  at  least  1  nesting  at¬ 
tempt  that  year). 

b.  The  maximum  carrying  capacity  (input  #8)  for  the  Fort  Stewart  RCW 
population  is  set  to  its  estimated  value  of  506  breeding  females  under  the 
most  recent  draft  version  of  the  Fort  Stewart  Endangered  Species  Man¬ 
agement  Plan. 

c.  The  pseudoextinction  threshold  (input  #9),  below  which  a  population  is 
considered  effectively  extinct,  was  set  to  5  breeding  females. 

d. The  target  population  size  (input  #10)  which  is  desirable  to  equal  or  ex¬ 
ceed  as  a  criterion  of  successful  management,  was  set  to  be  250  breeding 
females,  following  the  USFWS  guidelines  delineated  in  the  RCW  recovery 
plan. 

e.  The  maximum  possible  seasonal  fecundity  (input  #17)  was  set  to  the 
maximum  observed  size  of  4  nestlings  (of  both  sexes). 

f.  The  probability  that  a  fledgling  is  female  (input  #  18)  was  set  to  equal 
0.5,  following  the  observed  fledgling  sex  ratio  of  1:1. 
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g.  The  mean  fledgling  female  survival  rate,  and  its  associated  temporal 
and  sampling  standard  deviations  (inputs  #19,  20,  21). 

h.  The  mean  adult  female  survival  rate,  and  its  associated  temporal  and 
sampling  standard  deviations  (inputs  #22,  23,  24). 

i.  The  recurrence  rates,  and  the  effects  on  survival  and  fecundity,  of  hurri¬ 
cane  catastrophes  (inputs  #30  through  44)  are  somewhat  different  than 
those  in  the  previous  report,  although  their  derivations  are  basically  the 
same. 

2.  Parameters  whose  values  varied  among  disturbance  scenarios: 

a.  The  effective  sample  size  for  each  disturbance  scenario  (input  #3). 

b.  The  mean  proportion  of  successful  nests,  and  its  associated  temporal  and 
sampling  standard  deviations  (inputs  #11, 12, 13). 

c.  The  mean  number  of  fledglings  |  successful  nest,  and  its  associated  tem¬ 
poral  and  sampling  standard  deviations  (inputs  #14, 15  16). 

The  latter  parameters,  whose  values  varied  among  disturbance  scenarios,  are 
the  reproductive  statistics  previously  discussed  above. 

The  previous  study  using  PVAvES  (Melton  et  al.  2001)  showed  that  the  omission 
of  the  effects  of  hurricane  catastrophes,  in  a  population  viability  analysis  of  the 
Fort  Stewart  RCW  population,  results  in  a  considerable  underestimation  of  the 
likelihood  of  extinction.  Therefore,  hurricane  catastrophes  were  included  in  the 
present  analyses.  The  parameters  for  hurricane  occurrence  rates,  and  hurricane 
effects  on  survival  and  fecundity,  in  the  present  analysis  were  somewhat  differ¬ 
ent  than  those  used  in  the  previous  study.  The  parameter  input  values,  in  the 
present  study,  for  the  recurrence  rates  of  hurricanes  (inputs  #30,  33,  36,  39,  42) 
were  based  on  the  estimated  recurrence  rates  of  Saffir-Simpson  category  1,11, 
III,  IV,  and  V  strength  hurricanes  for  Fort  Stewart  (Hooper  and  McAdie  1995). 
The  parameter  values  for  percent  reduction  of  survival  rates,  and  population 
carrying  capacity,  by  Category  III  hurricanes  (inputs  #37,  38)  were  taken  from 
the  documented  effects  of  hurricane  Hugo  on  the  RCW  population  at  Francis 
Marion  National  Forest  in  North  Carolina  (Watson  et  al.  1995).  Hugo  was  a 
category  III  hurricane  when  it  hit  the  population.  In  the  absence  of  better  in¬ 
formation,  comparable  parameter  values  for  categories  I  ,  II,  IV,  and  V  hurri¬ 
canes  were  extrapolated  from  the  parameter  values  for  category  III  hurricanes. 
This  was  done  by  assuming  that  their  values  varied  in  direct  proportion  to  the 
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square  of  the  threshold  maximum  sustained  wind  speeds  for  each  category  as 
given  in  Hooper  and  McAdie  (1995),  i.e.,  varied  in  proportion  to  the  wind's  ki¬ 
netic  energy.  All  parameter  values  for  PVAvES  catastrophe  levels  6-10  (inputs  # 
45-59)  were  set  equal  to  0. 

Two  sets  of  PVAvES  outputs  were  reported  in  this  study,  for  each  disturbance 
scenario.  One  set  of  outputs  included  the  median  probabilities  of  pseudoextinc¬ 
tion  within  less  than  or  equal  to  10  years,  20  years,  and  100  years.  Medians, 
rather  than  means,  were  reported  because  the  underlying  distributions  were 
typically  highly  skewed  and/or  bimodal.  The  second  set  of  outputs  included  the 
probabilities  that  the  population  should  be  classed  as  VULNERABLE,  as 
ENDANGERED,  and  as  CRITICAL,  using  the  definitions: 

VULNERABLE:  The  probability  of  pseudoextinction  within  100  years  >  0.1 

ENDANGERED:  The  probability  of  pseudoextinction  within  20  years  >  0.2 

CRITICAL:  The  probability  of  pseudoextinction  within  10  years  >  0.5. 

These  definitions  are  taken  from  the  IUCN  extinction  risk  criteria,  as  proposed 
in  Mace  and  Stuart  (1994). 


Results 

Classification  of  Cluster  Sites  by  Military  Activity  Levels 

The  bivariate  distribution  of  clusters  based  on  observed  military  and  civilian  ac¬ 
tivity  is  presented  in  Figure  12.  This  activity  includes  all  civilian  and  military 
activity  as  described  in  Tables  2  through  5,  including  both  vehicle  and  dis¬ 
mounted  personnel  activity.  Most  clusters  did  not  show  exceedingly  high  levels 
of  either  kind  of  activity.  However,  three  clusters  (clusters  #7,  47,  and  57)  could 
be  readily  identified  as  having  unambiguously  high  relative  levels  of  military 
and  (to  a  lesser  degree)  civilian  activity.  These  three  clusters  were  designated  as 
"high”  activity  clusters,  and  the  rest  designated  the  "low"  activity  clusters.  This 
visually  patent,  but  somewhat  arbitrary,  classification  corresponds  to  a  post  hoc 
criterion  for  each  "high"  activity  cluster  site  of  mean  MILACT  >0.2,  and  mean 
CIVACT  >0.1. 
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Fledgling  Production  in  Each  Activity  Level 

Nest  productivity  statistics  for  each  activity  class  and  year,  along  with  their 
standard  errors  and  sample  sizes,  are  shown  in  Table  18.  The  3-year  means  of 
the  mean  proportion  of  successful  nests  (Table  19)  showed  no  significant  differ¬ 
ence  between  activity  classes  at  the  .05  level  (t  =  -0.6565,  df  =  4,  p  =  0.5473). 
The  3-year  means  of  the  mean  number  of  fledglings  |  successful  nest  (Table  19) 
did  differ  significantly,  at  the  5%  level,  between  activity  classes  (t  =-4.1587,  df  = 
4,  p=  0.0142). 


Mean  Proportion  of  Observation  Periods  with  Military 

Activity 


Figure  12.  Classification  of  Red-cockaded  Woodpecker  clusters  used  in  the  population  viability 
analysis,  based  on  2-  or  3-year  mean  levels  of  military  and  civilian  activity. 
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Table  18.  Nest  productivity  summary  statistics,  over  three  years,  for  Red-cockaded  Woodpecker 
clusters  in  “low”  and  “high”  disturbance  activity  classes. 


Nest  Productivity  ---  Descriptive  Statistics 

A)  Low  Activity: 

Year 

Proportion  Successful 

Mean  ±  Binomial  SE  (N)* 

Fledglings  |  Successful 

Mean  ±  SE  (N) 

1997 

0.9048  ±0.0641  (21) 

2.3158  ±0.1539  (19) 

1998 

0.761 9  ±0.0929  (21) 

2.1875  ±0.1638  (16) 

1999 

0.9047  ±0.0641  (21) 

2.2105  ±0.1806  (19) 

B)  High  Activity: 

Year 

Proportion  Successful 

Mean  ±  Binomial  SE  (N) 

Fledglings  |  Successful 

Mean  ±  SE 

1997 

0.6667  ±  0.2722  (3) 

1 .0000  ±  0.0000  (2) 

1998 

1 .0000  ±  0.0000  (3) 

1 .6667  ±  0.3333  (3) 

1999 

0.6667  ±0.2722  (3) 

1 .5000  ±  0.5000  (2) 

*“SE”  denotes  standard  error.  “N”  denotes  sample  size. 


Table  19.  Estimated  PVAvES  nest  productivity  parameters,  calculated  from  the  RCW  clusters  in 
“low”  and  “high”  disturbance  activity  classes. 


Military  Disturbance 
Category: 

Estimated  Nest  Productivity  Parameter 

Low 

High 

Mean  Proportion  Successful 

0.8571 

0.7778 

Temporal  SD  (standard  deviation) 

0.0345 

0.0000 

Sampling  SD 

0.0749 

0.2222 

Mean  Fledglings  |  Successful 

2.2379 

1 .3889 

Temporal  SD 

0.0000 

0.0000 

Sampling  SD 

0.1665 

0.3469 

Effective  Sample  Size 

63 

9 

Hypothetical  Disturbance  Scenarios 

Table  19  shows  the  3-year  means  of  the  means  in  Table  18,  along  with  their 
temporal  (over  3  years)  and  sampling  error  standard  deviations,  and  associated 
effective  sample  sizes.  Table  20  presents  the  nest  productivity  parameter  values, 
and  effective  sample  sizes,  for  all  9  hypothetical  disturbance  scenarios  used  as 
inputs  to  the  PVAvES  program.  Scenario  values  were  calculated  as  weighted 
averages  of  the  values  for  "low”  and  "high”  disturbance  categories  shown  in  Ta¬ 
ble  19. 
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Table  20.  Parameterization  of  nine  hypothetical  disturbance  scenarios,  derived  as  weighted 
averages  of  the  parameter  estimates  in  Table  19. 


Militarv  Disturbance  Scenarios: 

Hypothetical  proportion  of  RCW  habitat  in  High  Military  Disturbance 
category 

Nest  Productivity 

Parameter 

PVAvES 

Input  ID 

0 

0.125 

0.250 

0.375 

0.500 

0.625 

0.750 

0.875 

1 

Mean  Proportion  Suc¬ 
cessful 

#  11 

0.8571 

0.8472 

0.8373 

0.8274 

0.8175 

0.8075 

0.7976 

0.7877 

0.7778 

Temporal  SD  (standard 
deviation) 

#12 

0.0345 

0.0302 

0.0258 

0.0215 

0.0172 

0.0129 

0.0086 

0.0043 

0.0000 

Sampling  SD 

#13 

0.0749 

0.0933 

0.1118 

0.1302 

0.1486 

0.1670 

0.1854 

0.2038 

0.2222 

Mean  Fledglings  | 

Successful 

#14 

2.2379 

2.1318 

2.0257 

1.9195 

1.8134 

1.7073 

1.6012 

1 .4950 

1 .3889 

Temporal  SD 

#  15 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Sampling  SD 

#16 

0.1665 

0.1890 

0.2116 

0.2342 

0.2567 

0.2793 

0.3018 

0.3244 

0.3469 

Effective  Sample  Size 

#3 

63.00 

56.25 

49.50 

42.75 

36.00 

29.25 

22.50 

15.75 

9.00 

Population  Viability  Analysis  using  PVAvES  1.0 

Figure  13  illustrates  the  relationships  of  the  median  probabilities  of  pseudoex¬ 
tinction  within  100  years,  20  years,  and  10  years,  to  the  proportion  of  RCW  habi¬ 
tat  in  the  "high  activity"  category.  The  median  probability  of  pseudoextinction 
within  100  years  showed  a  sigmoid  relationship  to  the  proportion  of  "high"  activ¬ 
ity  habitat  (Figure  13A),  increasing  steeply  for  proportions  greater  than  25  per¬ 
cent  A  probability  of  pseudoextinction,  within  100  years,  of  100  percent  could  not 
be  ruled  out  (based  on  95  percent  confidence  intervals)  for  proportions  of  "high" 
activity  habitat  greater  than  50  percent  The  median  probability  of  pseudoextinc¬ 
tion  within  20  years  (Figure  13B)  was  generally  low,  but  increased  steeply  when 
the  proportion  of  "high"  activity  habitat  was  greater  than  75  percent.  The  me¬ 
dian  probability  of  pseudoextinction  within  10  years  (Figure  13C)  was  very  low 
for  all  disturbance  scenarios. 

Figure  14  illustrates  the  relationships  of  the  probabilities  that  the  Fort  Stewart 
RCW  population  should  be  classified  by  the  risk  categories  VULNERABLE, 
ENDANGERED,  and  CRITICAL,  to  the  proportion  of  RCW  habitat  in  the  "high 
activity"  category.  The  probability  that  the  class  =VULNERABLE  (Figure  14A) 
was  significantly  greater  than  0  (based  on  95  percent  confidence  intervals)  even 
when  the  proportion  of  "high"  activity  habitat  was  0.  It  showed  a  concave  in¬ 
crease  with  increasing  proportions  of  "high"  activity  habitat,  and  probabilities  of 
90  percent  or  more  could  not  be  ruled  out  for  proportions  of  "high"  activity  habi¬ 
tat  greater  than  75  percent.  The  probability  that  the  class  =  ENDANGERED 
(Figure  14B)  was  low  when  the  proportion  of  "high"  activity  habitat  was  0,  and 
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showed  a  convex  increase  with  increasing  proportions  of  "high"  activity  habitat. 
It  was  significantly  greater  than  0  for  proportions  of  "high"  activity  habitat 
greater  than  or  equal  to  50  percent.  The  probability  that  the  class  =  CRITICAL 
(Figure  14C)  was  very  low  for  all  disturbance  scenarios. 
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A 


— ♦ — Bootstrap  Estimate 
a  Lower  95%CI 
x  Upper  95%CI 


B 


-♦ — Bootstrap  Estimate 
a  Lower95%CI 
x  Upper  95%CI 


c 


-♦—Bootstrap  Estimate 
A  Lower95%CI 
x  Upper  95%CI 


Proportion  "High  Activity"  Habitat 


A,  within  100  years;  B,  within  20  years;  and  C,  within  10  years,  when  different  hypothetical  proportions 
of  the  habitat  are  assumed  to  be  “high  activity”  habitat.  Error  bars  denote  1  bootstrap  standard  error. 
“Cl”  denotes  bootstrap  percentile  confidence  interval. 


Figure  13.  Estimated  pseudoextinction  probabilities  for  the  Red-cockaded  Woodpecker 
population  on  Fort  Stewart. 
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Proportion  "High  Activity"  Habitat 
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Proportion  "High  Activity"  Habitat 


0  0.2  0.4  0.6  0.8  1 

Proportion  "High  Activity"  Habitat 


A,  VULNERABLE;  B,  ENDANGERED;  and  C,  CRITICAL  when  different  hypothetical  proportions  of  the 
habitat  are  assumed  to  be  “high  activity”  habitat.  Error  bars  denote  1  bootstrap  standard  error.  “Cl”  de¬ 
notes  bootstrap  percentile  confidence  interval. 

Figure  14.  Estimated  probabilities  that  the  Red-cockaded  Woodpecker  population  on  Fort 
Stewart  should  be  classified  as  vulnerable,  endangered,  and  critical. 
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Discussion 

Division  of  the  RCW  cluster  data  set  into  "high”  and  "low”  activity  classes  proved 
to  be  a  useful  device  for  the  present  population  viability  analysis,  allowing  crea¬ 
tion  of  the  nine  hypothetical  disturbance  activity  scenarios.  Only  3  clusters  out 
of  the  24  shown  in  Figure  12  (and,  in  fact,  only  3  out  of  all  51  clusters  with  2  or 
more  years  with  military-civilian  activity  observations)  stood  out  as  "high”  activ¬ 
ity  clusters  —  a  very  small  sample  size.  However,  the  fact  that  nest  success,  and 
fledglings  |  success,  were  lower  in  the  "high”  than  in  the  "low”  activity  clusters 
(significantly  so  for  fledglings  |  success)  suggests  that  some  correlate  of  military 
or  civilian  activity  may  indeed  reduce  nest  productivity  in  these  "high"  activity 
clusters.  Thus  the  question  "What  would  happen  in  terms  of  risk  to  the  popula¬ 
tion  if  more  RCW  habitat  on  the  installation  were  subject  to  "high”  activity  lev¬ 
els?”  is  a  logical  one  to  ask. 

When  interpreting  the  results  presented  here,  it  is  best  to  keep  in  mind  that  the 
estimates  reported  for  the  population  viability  parameters  are  not  necessarily 
conservative  estimates.  Although  demographic,  environmental,  catastrophic, 
and  sampling  variation  in  life-history  parameters  are  incorporated  into  the 
PVAvES  model,  no  model  is  likely  to  include  all  potential  sources  of  variation, 
and  any  such  unaccounted  sources  would  most  likely  tend  to  further  reduce 
population  viability  estimates.  Also,  survival  rates  of  breeding  females  are,  in 
the  present  study,  assumed  to  be  the  same  for  "low”  and  "high"  activity  habitats. 
Therefore,  the  reported  parameter  estimates  should  not  be  taken  too  literally. 
Comparisons  among  the  hypothetical  disturbance  scenarios,  however,  are  more 
likely  to  be  valid. 

With  this  caveat  in  mind,  the  present  results  suggest  that  if  a  substantial  pro¬ 
portion  (e.g.,  over  50  percent  of  RCW  habitat  on  the  installation  were  subject  to 
"high”  military-civilian  activity,  this  would  have  a  significant  negative  effect  on 
the  viability  of  the  Fort  Stewart  RCW  population.  Such  a  condition  would  in¬ 
crease  the  population's  probability  of  extinction  (especially  in  the  100-year  long 
term),  and  concomitantly  increase  the  likelihood  that  the  population  should  be 
classified  as  VULNERABLE  and  also  as  ENDANGERED. 

The  present  proportion  of  RCW  habitat  on  Fort  Stewart  that  is  subject  to  "high” 
military-civilian  activity  is  not  known  with  certainty.  Out  of  the  sample  of  51 
cluster  sites  having  at  least  2  years  of  observation  data  on  military  and  civilian 
activity,  only  3  clusters  (the  "high"  activity  clusters  shown  in  Figure  12),  or  6 
percent  of  the  total  sample,  had  mean  values  of  MILACT  >  0.2.  If  we  take  a 
value  of  less  than  10  percent  to  roughly  represent  the  proportion  of  "high” 
activity  RCW  breeding  habitat  on  Fort  Stewart,  then  the  results  suggest  that  the 
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viability  of  the  RCW  population  is,  at  present,  essentially  unaffected  by  "high" 
levels  of  military  activities.  Whether  this  situation  will  continue,  when  training 
restrictions  are  revised  in  accordance  with  the  1996  RCW  management 
guidelines,  needs  to  be  monitored. 

Even  if  all  of  the  RCW  breeding  habitat  on  the  installation  were  in  the  "low”  ac¬ 
tivity  class,  Figure  14A  shows  that  the  Fort  Stewart  RCW  population  would  have 
a  24  to  48  percent  chance  that  its  viability  status  should  be  considered 
VULNERABLE,  given  the  model's  assumptions  and  the  present  input  data.  The 
population,  while  not  presently  in  extremis,  should  not  be  considered  safe  from 
extinction  within  the  next  100  years,  even  under  a  low  military-civilian  distur¬ 
bance  regime.  Continued  monitoring,  and  management,  of  the  Fort  Stewart 
RCW  population  appears  to  be  warranted. 
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8  Conclusions  and  Recommendations 


Fort  Stewart  Results 

Results  from  this  study  on  Fort  Stewart,  Georgia  during  1997-1999  indicate  that 
demographic  factors  (e.g.,  group  size  and  prior  reproductive  success)  have  more 
effect  on  RCW  reproductive  success  than  habitat  and/or  disturbance  from  human 
activities.  Observations  of  human  activity  at  RCW  sites  suggest  that  the  prob¬ 
ability  of  disturbance  from  military  training  activities  in  clusters  is  relatively 
low  in  the  majority  of  RCW  cluster  sites  on  Fort  Stewart.  Flowever,  data  from  a 
small  number  of  clusters  in  high-traffic  areas  on  the  installation  suggest  that 
disturbance  exceeding  certain  levels  of  activity  could  be  detrimental  to  RCW  re¬ 
productive  success.  Population  viability  modeling  indicates  that  at  the  present 
time  potential  disturbance  effects  in  this  small  proportion  of  the  population  have 
negligible  effect  on  the  viability  of  the  Fort  Stewart  RCW  population. 

These  results  indicate  that  current  Fort  Stewart  management  practices  are  suc¬ 
cessfully  mitigating  variance  in  reproductive  parameters  attributable  to  effects 
of  habitat.  Several  factors  in  previously  published  studies  (Carter,  Walters,  and 
Doerr  1995;  DeLotette,  Epting,  and  DeMuth  1995)  have  related  cluster  occu¬ 
pancy  and  reproductive  success  to  numbers  of  available  cavities  and  midstory 
encroachment.  This  study  did  not  find  any  significant  association  between  these 
habitat  factors  and  cluster  occupancy  or  reproductive  success  in  monitored  clus¬ 
ters.  These  findings  are  not  surprising  given  the  aggressive  habitat  manage¬ 
ment  in  practice  at  Fort  Stewart  over  the  past  several  years.  These  results  indi¬ 
cate  that  aggressive  management  can  minimize  habitat  as  a  limiting  factor  on 
RCW  populations.  Continuing  these  management  practices  in  all  primary  and 
supplemental  clusters  will  continue  to  reduce  potential  limits  on  RCW  popula¬ 
tions  due  to  habitat  features. 

Demographic  factors  (group  size  and  prior  reproductive  success)  have  the  most 
discernible  relationship  to  RCW  reproductive  success  on  Fort  Stewart.  Demo¬ 
graphic  limitations  on  populations  are  more  difficult  to  ameliorate  through  direct 
management  intervention  than  habitat  management  practices.  Flowever,  pro¬ 
viding  the  necessary  conditions  such  as  adequate  cavity  availability  and  mini¬ 
mizing  cluster  isolation  will  be  supportive  of  desirable  demographic  profiles  in 
the  Fort  Stewart  population. 
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During  the  time  period  of  this  study  (1997-1999)  training  restrictions  in  effect  on 
Fort  Stewart  were  in  accordance  with  the  1994  "Management  Guidelines  for 
RCWs  on  Army  Installations.”  These  restrictions  essentially  prohibited  mecha¬ 
nized  maneuver  training  activities  within  200  feet  of  cavity  or  start  trees  except 
on  maintained  roads  and  trails.  Adherence  to  these  restrictions  was  reflected  in 
observations  of  training  activity  in  proximity  to  RCW  clusters  during  1997-1999. 
The  minimal  effect  of  maneuver  training  activities  during  this  period  was  ob¬ 
served  in  the  context  of  operations  under  these  training  restrictions. 

Prior  to  the  2001  RCW  breeding  season,  Fort  Stewart  will  implement  its  instal¬ 
lation  Endangered  Species  Management  Plan  (ESMP)  in  accordance  with  the 
revised  1996  "Management  Guidelines  for  RCWs  on  Army  Installations."  Re¬ 
strictions  on  maneuver  training  operations  in  proximity  to  RCW  clusters  will  be 
significantly  reduced  under  the  implemented  ESMP.  Military  vehicle  traffic  will 
be  unrestricted  within  cluster  sites  except  within  a  50-ft  buffer  around  cavity 
and  start  trees. 

This  study  provides  the  baseline  data  necessary  to  evaluate  potential  effects  of 
implementing  revised  training  restrictions  on  RCW  populations  on  Fort  Stewart, 
Georgia.  This  study  identified  monitored  clusters  that  consistently  receive 
higher  levels  of  military  training.  These  likely  will  be  the  clusters  that  are  at 
higher  risk  from  potential  effects  of  training  under  the  reduced  restrictions. 
Monitoring  data  collected  by  installation  biologists  in  clusters  with  relatively 
high  levels  of  training  activity  should  be  assessed  annually  in  comparison  with 
baseline  data  from  the  clusters  during  1997-1999  and  relative  to  clusters  identi¬ 
fied  at  low  risk  from  military  activity  within  the  monitoring  year.  This  compari¬ 
son  should  be  valid  unless  the  Fort  Stewart  training  mission  significantly 
changes  or  distribution  of  training  activity  significantly  changes. 

Based  on  results  of  this  study  it  is  not  recommended  that  intensive  site-specific 
monitoring  of  training  activity  be  continued.  The  bases  for  this  recommendation 
are  that  (1)  the  proportion  of  observed  training  activity  is  relatively  low,  which  is 
problematic  for  establishing  statistical  relationships  between  frequency  of  train¬ 
ing  and  measures  of  RCW  population  parameters,  and  (2)  as  noted  above,  rela¬ 
tive  risk  of  exposure  to  training  activity  has  been  identified  for  monitored  clus¬ 
ters  in  this  study,  assuming  no  significant  change  in  the  conduct  of  training 
operations  on  Fort  Stewart.  Use  of  RF MSS  scheduling  data  is  recommended  to 
assess  the  overall  level  and  distribution  of  training  activity  during  the  baseline 
years  1998  and  1999  and  future  years  of  population  monitoring. 

Data  presented  in  this  report  show  a  weak  association  with  observed  training  at 
the  site  level,  which  is  likely  an  artifact  of  the  overall  low  frequency  of  observed 
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activity  at  the  site  level  and  the  low  "resolution"  of  the  RFMSS  scheduling  vari¬ 
ables  used  in  analyses  in  this  report.  However,  comparison  of  RFMSS  data 
among  years  should  provide  an  index  of  change  in  conduct  of  training  operations 
on  the  installation. 

It  is  recommended  that  site-specific  monitoring  of  training  activity  should  be  im¬ 
plemented  only  if  (1)  significant  negative  changes  in  monitored  population  pa¬ 
rameters  are  observed  between  data  obtained  during  baseline  years  1997-1999 
and  data  in  later  years,  or  (2)  the  training  mission,  training  intensity,  or  conduct 
of  training  operations  on  Fort  Stewart  significantly  changes  relative  to  scheduled 
activity  during  1997-1999. 


Use  in  Endangered  Species  Risk  Assessments 

Results  of  this  study  cannot  be  extrapolated  directly  to  RCW  populations  on 
other  military  installations.  Results  presented  here  are  largely  dependent  on  the 
combination  of  habitat,  demographic,  and  training  parameters  that  are  unique  to 
Fort  Stewart.  However,  the  approaches  taken  in  this  study  to  evaluate  risk  to 
RCW  populations  due  to  training  activities  are  applicable  to  other  installations 
with  RCW  populations  and  risk  assessment  for  other  endangered  species  popula¬ 
tions  on  military  installations. 

The  key  elements  of  this  approach  are: 

2.  Evaluation  of  exposure  probability  to  risk  factors  of  concern. 

3.  Evaluation  of  risk  factors  in  context  of  all  potential  limiting  factors  on  the  popu¬ 
lation  of  concern. 

4.  Evaluation  of  effects  of  potential  risk  factors  at  the  population  level. 

The  importance  of  the  first  element  is  illustrated  from  data  obtained  on  Fort 
Stewart  and  in  many  ways  is  the  most  difficult  to  evaluate  in  the  context  of 
training  activities  on  military  installations.  The  data  for  Fort  Stewart  clearly 
show  that  the  risk  of  exposure  to  maneuver  training  activities  on  Fort  Stewart  is 
relatively  low  for  the  majority  of  RCW  clusters.  Experimental  approaches  can 
evaluate  levels  of  effect  on  populations,  but  cannot  determine  potential  exposure 
to  any  particular  level  of  effect.  Data  for  Fort  Stewart  suggest  a  level  of  signifi¬ 
cant  effect.  However,  the  majority  of  the  population  is  not  exposed  to  this  level 
of  effect.  I  n  the  context  of  military  operations,  determining  exposure  risk  is  diffi¬ 
cult  given  the  dynamic  and  complex  spatial  and  temporal  characteristics  of  ma¬ 
neuver  training  activities.  In  this  study,  exposure  to  training  activity  was 
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evaluated  by  sampling  training  activity  at  cluster  sites.  This  approach  is  appli- 
cablefor  any  risk  assessment  to  determine  probability  of  site-specific  exposure  to 
training  activity.  However,  in  large-scale  applications,  this  approach  is  expen¬ 
sive,  time  consuming,  and  personnel  intensive.  I  mplementing  this  approach  un¬ 
der  all  circumstances  may  allocate  limited  management  resources  away  from  ac¬ 
tivities  that  directly  benefit  populations  of  concern.  In  many  cases  exposure  to 
risk  factors  from  maneuver  training  activity  may  be  most  efficiently  determined 
by  a  relative  classification  of  risk  based  on  expert  opinion  and  readily  available 
data  from  sources  such  as  RFMSS  scheduling  data.  Data  from  assigned  risk  lev¬ 
els  can  then  be  assessed  in  relation  to  other  factors  known  or  predicted  to  limit 
populations. 

In  RCWs,  several  factors  have  been  identified  that  affect  reproduction  and  clus¬ 
ter  occupancy  including  both  demographic  and  habitat  parameters.  This  report 
analyzes  potential  training  effects  in  relation  to  these  other  parameters;  in  fact, 
demographic  factors  were  the  best  predictors  of  some  aspects  of  RCW  reproduc¬ 
tive  success  and  cluster  occupancy.  This  result  emphasizes  the  importance  of 
installation  managers  being  aware  of  the  biological  factors  potentially  limiting 
populations  under  their  management  concern  and  of  tailoring  monitoring  plans 
to  acquire  the  relevant  data  necessary  to  evaluate  all  potential  sources  of  popula¬ 
tion  limits.  For  installations  implementing  the  Army  1996  "Management  Guide¬ 
lines  for  the  Red-cockaded  Woodpecker  on  Army  I  nstallations"  estimates  of  criti¬ 
cal  demography  parameters  will  be  available.  However,  based  on  data 
availability  on  Fort  Stewart,  special  emphasis  should  be  placed  on  follow-up 
monitoring  of  banded  populations.  Emphasis  on  intensive  habitat  data  collection 
is  probably  not  warranted  if  installations  are  adequately  implementing  pre¬ 
scribed  burning  and  cavity  augmentation  programs.  Data  for  Fort  Stewart  sug¬ 
gest  that  habitat  is  not  an  important  factor  in  clusters  with  adequate  midstory 
control  and  cavity  availability. 

The  most  important  aspect  of  this  study's  research  approach  was  to  evaluate 
effects  of  observed  population  variability  at  the  population  level.  Although 
effects  at  the  individual  level  have  consequences  in  a  regulatory  and  compliance 
context,  ultimately  it  is  the  effect  at  the  population  level  that  is  significant  to 
recovery  objectives  and  prescriptions  for  protection  requirements.  The 
population  viability  approach  presented  in  this  report  provides  a  probabilistic 
evaluation  of  extinction  risk  and  population  recovery  potential.  This  is 
particularly  relevant  for  installations  establishing  population  objectives  for 
installation  ESMPs.  This  approach  allows  managers  to  evaluate  the  effect  of 
observed  demographic  variability  at  the  population  level.  Although  evaluating 
installation  population  monitoring  data  in  viability  models  cannot  account  for 
the  cause  of  change  in  extinction  risk  or  recovery  objectives,  significant  changes 
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in  extinction  risk  or  recovery  objectives  can  be  the  trigger  for  implementing  more 
intensive  monitoring  or  research  programs  to  determine  causal  associations. 
Based  on  results  presented  in  this  report  it  is  recommended  that  installations 
evaluate  population  extinction  risk  and  recovery  probabilities  for  populations  of 
management  concern  and  formally  identify  change  thresholds  for  implementing 
more  intensive  programs  to  identify  and  mitigate  risk  factors. 
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Appendix:  PVAvES  Input  Data  File 

Printout  of  the  PVAvES  formatted  input  data  file,  used  in  this  study  to  param¬ 
eterize  all  nine  hypothetical  disturbance  scenarios.  Parameter  values  denoted  as 
"VARIED"  differed  among  the  scenarios,  as  shown  in  Table  20.  All  other  pa¬ 
rameter  values  were  exactly  the  same  for  all  scenarios. 


************  INPUT  DATA  FOR  PROGRAM  PVAvES  ************ 


***  INITIALIZATION  PARAMETERS  *** 


1)  Seedl  =  123456 

2)  Seed2  =  765432 

3)  Effective  Sample  Size  from  Demographic  Data  -  VARIED 

4)  Number  of  Bootstrap  Iterations  -  2000 

5)  Size  of  Approximated  Sampling  Distribution  -  10000 

6)  Coverage  of  Confidence  Intervals  =  0.95 


***  POPULATION  AND  THRESHOLD  PARAMETERS  *** 


7)  Initial  population  size  -  165 

(breeding  females) 

8)  Maximum  population  carrying  capacity  =  506 

(breeding  females,  in  the  absence  of  catastrophes) 
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9)  Extinction  threshold  - 

10)  Target  population  size  after  100  years  = 


***  DEMOGRAPHIC  PARAMETERS  *** 


5 

250 


11)  Mean  proportion  of  successful  nests  = 

12)  with  temporal  standard  deviation  = 

13)  and  sampling  standard  deviation  = 

14)  Mean  seasonal  fecundity  of  successful  nests  (fledglings  of  both  sexes) 

15)  with  temporal  standard  deviation  = 

16)  and  sampling  standard  deviation  = 

17)  Maximum  possible  seasonal  fecundity  = 

(fledglings  of  both  sexes) 

18)  Probability  that  a  fledgling  is  female  = 

(based  on  fledgling  sex  ratio) 

19)  Mean  fledgling  female  survival  rate  (to  next  breeding  season)  - 

20)  with  temporal  standard  deviation  = 

21)  and  sampling  standard  deviation  - 

22)  Mean  adult  female  annual  survival  rate  = 

23)  with  temporal  standard  deviation  = 

24)  and  sampling  standard  deviation  = 


VARIED 

VARIED 

VARIED 

VARIED 

VARIED 

VARIED 

4 

0.5 

0.3901 

0.0000 

0.0624 

0.7092 

0.0402 

0.0572 


***  IMMIGRATION  PARAMETERS  *** 

itititititititititit'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

25)  Immigration  constant  or  stochastic  (poisson) ?  =  0 

(Enter  0  for  CONSTANT,  1  for  STOCHASTIC  ) 

26)  Mean  immigration  rate  (breeding  females  per  year)  -  0 

(must  be  an  integer  for  CONSTANT  immigration) 
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********  CATASTROPHE  PARAMETER  VALUES  ******** 


27)  Catastrophes  are  included  in  model,  and  reduce  annual  survival  rates?  -  1 

(Enter  0  or  1 :  0  for  NO,  1  for  YES) 

28)  Catastrophes  also  reduce  carrying  capacity?  -  1 

(Enter  0  or  1 :  0  for  NO,  1  for  YES) 

29)  Regeneration  time  of  carrying  capacity  (years)  =  100 

(Time  for  carrying  capacity  to  grow  from  1%  to  99%  of  its  maximum  value) 

Category  1  Catastrophes: 

30)  Rate  of  occurrence  (annual)  -  0.04 

31)  Proportional  reduction  in  survival  rates  =  0.31 

32)  Proportional  reduction  in  carrying  capacity  =  0.26 

Category  2  Catastrophes: 

33)  Rate  of  occurrence  (annual)  =  0.009090 

34)  Proportional  reduction  in  survival  rates  =  0.52 

35)  Proportional  reduction  in  carrying  capacity  -  0.44 

Category  3  Catastrophes: 

36)  Rate  of  occurrence  (annual)  =  0.00357 

37)  Proportional  reduction  in  survival  rates  =  0.70 

38)  Proportional  reduction  in  carrying  capacity  -  0.59 

Category  4  Catastrophes: 

39)  Rate  of  occurrence  (annual)  -  0.00095 

40)  Proportional  reduction  in  survival  rates  =  0.97 

41)  Proportional  reduction  in  carrying  capacity  =  0.82 

Category  5  Catastrophes: 

42)  Rate  of  occurrence  (annual)  -  0.000018 

43)  Proportional  reduction  in  survival  rates  -  0.99 

44)  Proportional  reduction  in  carrying  capacity  =  0.99 

Category  6  Catastrophes: 
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45)  Rate  of  occurrence  (annual)  -  0.0000 

46)  Proportional  reduction  in  survival  rates  =  0.0000 

47)  Proportional  reduction  in  carrying  capacity  =  0.0000 

Category  7  Catastrophes: 

48)  Rate  of  occurrence  (annual)  -  0.0000 

49)  Proportional  reduction  in  survival  rates  =  0.0000 

50)  Proportional  reduction  in  carrying  capacity  =  0.0000 

Category  8  Catastrophes: 

51)  Rate  of  occurrence  (annual)  =  0.0000 

52)  Proportional  reduction  in  survival  rates  =  0.0000 

53)  Proportional  reduction  in  carrying  capacity  =  0.0000 

Category  9  Catastrophes: 

54)  Rate  of  occurrence  (annual)  =  0.0000 

55)  Proportional  reduction  in  survival  rates  =  0.0000 

56)  Proportional  reduction  in  carrying  capacity  -  0.0000 

Category  10  Catastrophes: 

57)  Rate  of  occurrence  (annual)  -  0.0000 

58)  Proportional  reduction  in  survival  rates  =  0.0000 


59) 


Proportional  reduction  in  carrying  capacity  =  0.0000 
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